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Abstract 

The  NAVSTAR  Global  Positioning  System  (GPS)  is  currently  used  in  many  appli¬ 
cations  requiring  precise  positioning  data.  Improving  the  precise  positioning  information 
requires  the  removal  of  errors  that  perturb  the  received  signals.  The  errors  introduced  by 
multiple  propagation  channels,  termed  multipath,  axe  not  easily  removed.  These  channels 
are  caused  by  reflective  surfaces  near  the  receiver.  As  such,  multipath  is  uncorrelated 
between  receivers  and,  thus,  cannot  be  removed  through  differencing  techniques. 

This  thesis  investigates  a  GPS  code  tracking  loop  design  which  uses  maximum  like¬ 
lihood  (ML)  estimation  to  determine  amplitude  and  phase  information  of  the  multipath 
signal  which  are  used  to  adjust  code  tracking  to  account  for  multipath  effects.  Analysis  of 
the  operations  that  govern  this  design  for  the  case  of  a  single  reflection  shows  that  it  has 
no  steady  state  tracking  error.  Results  of  simulations  indicate  that  the  code  tracking  loop, 
in  conjunction  with  the  MLE,  mitigate  the  effects  of  multipath  emd  improves  code  track¬ 
ing  performance  over  the  narrow  correlator  NCDLL  for  most  scenarios  analyzed.  Overall 
results  of  simulations  indicate  that  the  implementation  of  the  maximum  likelihood  esti¬ 
mator  (MLE)  in  conjimction  with  the  code  tracking  loop  has  the  potential  to  enhance 
code  tracking  performance  over  that  offered  by  the  narrow  correlator  NCDLL  in  a  GPS 
environment. 
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A  GPS  CODE  TRACKING  RECEIVER  DESIGN 


FOR  MULTIPATH  MITIGATION 
USING  MAXIMUM  LIKELIHOOD  ESTIMATION 

I.  Introduction 

1.1  Overview 

The  NAVSTAR  Global  Positioning  System  (GPS)  is  currently  used  in  many  applica¬ 
tions  in  both  the  military  and  civilian  community  to  obtain  precise  positioning  data.  The 
popularity  of  GPS  continues  to  grow,  in  part,  because  of  its  eiccessibility  to  both  commu¬ 
nities.  Because  of  its  wide  use  and  acceptance,  applications  requiring  increased  precision 
continue  to  be  discovered. 

Improving  the  collection  of  precise  ranging  and  positioning  information  requires  the 
removal  of  errors  that  perturb  the  received  signals.  Many  of  the  errors  may  be  reduced 
or  completely  removed  by  differencing,  a  common  technique  used  to  remove  measurement 
errors  that  are  closely  correlated  between  two  or  more  neeurby  receivers.  Differencing, 
however,  does  not  reduce  errors  introduced  by  multiple  propagation  channels  commonly 
known  as  multipath,  because  the  multipath  signals  2ire  unique  for  each  receiver  and  thus 
uncorrelated. 

The  errors  introduced  as  a  result  of  multipath  can  be  significzmt.  As  such,  a  con¬ 
siderable  amount  of  research  has  been  directed  towards  the  mitigation  of  multipath.  This 
mitigation,  however,  is  complicated  by  the  feict  that  the  presence  of  multipath  is  not  read¬ 
ily  known.  In  addition,  the  effects  of  multipath  vary  with  time,  particularly  in  a  dynamic 
environment  such  as  that  encountered  by  an  aircraft  in  flight. 

1.2  Problem 

Multipath  or  multiple  propagation  paths  occur  when  the  transmitted  GPS  signal 
reflects  off  smrounding  objects,  such  as  the  ground,  buildings,  the  fuselage  or  wings  of 
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GPS  Satellite 


Figure  1  Multipath  signal  reflections 


an  aircraft,  emd  mountains,  prior  to  reception  at  the  receiver  as  illustrated  in  Figure  1. 
The  reflected  signals  experience  a  larger  propagation  delay  relative  to  the  direct  path  or 
line-of-sight  (LOS)  signal  because  of  the  additional  distance  traveled.  Under  certeun  con¬ 
ditions,  these  reflections  induce  errors  in  the  code  or  phase  synchronization  process  of  the 
receiver  that  translate  in  tiim  from  ranging  errors  into  positioning  errors.  Because  of  the 
correlation  characteristics  of  the  spreading  code,  for  traditional  code  tracking  loops,  mul¬ 
tipath  signals  with  delays  greater  than  1.5Tc,  where  Tc  is  the  chip  period  of  the  spreading 
code,  have  no  impact  on  the  pseudo-range  measurement  [19].  This  delay  interval  equates 
to  signals  that  axe  reflected  within  450  meters  of  the  receiver.  However,  multipath  signals 
within  450  meters  (delays  less  than  1.5rc)  can  introduce  significant  errors  into  the  ranging 
measurements.  This  thesis  investigates  a  GPS  code  trax:king  receiver  which  uses  maximum 
likelihood  (ML)  estimation  to  mitigate  the  effects  of  multipath. 
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1.3  Summary  of  Current  Knowledge 

To  investigate  the  performance  of  code  tracking  loops  in  the  presence  of  multipath, 
one  must  understand  the  composition  of  the  transmitted  signal  and  the  multipath  compo¬ 
sition  of  the  received  signal.  The  following  section  discusses  these  issues  and  also  describes 
the  correlation  operation  required  of  GPS  receivers. 

1.3.1  The  GPS  Signal  GPS  signals  are  composed  of  three  basic  elements:  the 
wide  band  spreading  code,  the  carrier  signal,  and  the  navigation  message  or  navigation 
code.  The  wide  band  spreading  code  possesses  pseudo  random  properties  which  essentially 
appear  as  noise  within  the  signal  bandwidth.  Upon  successful  correlation  with  a  locally 
generated  replica,  the  spreading  code  provides  a  processing  gain  which  eillows  signals  to 
be  transmitted  near  or  below  the  noise  floor.^  GPS  incorporates  two  types  of  spreading 
codes:  course/acquisition  (C/A)  code  and  precision  code  (P-code).  C/A  code  is  accessible 
to  the  civilian  community  while  P-code  is  intended  for  military  use  only.  The  C/A  code 
has  a  chip  rate  or  bit  rate  of  1/Tc  =  1.023  MHz  and  a  code  period,  N  =  1023  chips; 
while  the  P-code  has  a  chip  rate,  l/Tc  =  10.23  MHz  and  a  code  period,  N  one  week. 
The  null  to  null  bandwidth  of  eaxrh  code  is  2.046  MHz  and  20.46  MHz,  respectively.  The 
GPS  satellite  transmits  right  hand  circularly  polarized  (RHCP)  carrier  signals  within  the 
L-band  at  two  carrier  frequencies,  LI  =  1575.42  MHz  and  L2  =  1227.6  MHz.  The  C/A 
code  and  P-code  are  both  mixed  with  the  LI  carrier  while  only  the  P-code  is  transmitted 
via  the  L2  carrier.  Prom  this  point  forward,  it  is  assumed  that  only  C/A  code  is  being 
tracked  unless  otherwise  stated. 

The  navigation  message  provides  satellite  ephemerides,  system  time,  correction  data, 
satellite  clock  behavior  data,  and  status  messages,  necessary  to  ensure  correct  positioning 
solutions.  This  data  message  has  a  bit  rate  of  50  hits /sec.  The  navigation  data  is  combined 
with  the  pseudo  random  noise  (PRN)  ranging  code  or  spreading  code  by  Modulo-2  addition 
and  can  be  represented  as  a  string  of  ones  and  zeros.  This  bit  stream  is  then  modulated 
onto  the  carrier  using  binary  phase  shift  keying  (BPSK),  which  results  in  instemtaneous 

^The  processing  gain,  which  is  the  ratio  between  the  spreading  code  chip  rate  eind  the  data  modulation 
bit  rate,  for  GPS  C/A  code  is  2e4  =:  43(jB 
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phase  changes  of  the  carrier  by  180  degrees.  This  bandwidth  spreading  process  is  known  as 
direct-sequence  spread  spectrum  (DS/SS).  The  null-to-nuU  bandwidth  of  the  transmitted 
signal  is  2.046  MHz,  the  bandwidth  of  the  spreading  code,  and  is  centered  about  the  LI 
carrier  frequency.  The  transmitted  signal  can  be  represented  as 

s(t)  =  (t)  c  (t)  cos(27r/ot  +  0)  (1) 

where  P  is  the  power  of  the  signal,  c{t)  is  the  spreading  code,  m{t)  is  the  navigation 
message,  /o  is  the  carrier  frequency  in  Hz,  and  6  is  the  transmitted  signal  phase. 

1,S.2  Multipath  Propagation.  The  influence  of  the  reflected  signals  depends 
largely  on  their  amplitudes  and  delay  relative  to  the  direct  path  signed.  The  direct  signed 
will  have  more  power  than  the  reflected  signals,  unless  there  are  obstructions  in  the  LOS 
that  reduce  the  power  level  of  the  direct  signal.  This  is  attributed,  in  part,  to  the  longer 
distance  that  the  reflected  signals  must  travel,  which  causes  additional  free  space  loss.  More 
significantly,  there  is  typically  attenuation  caused  by  the  reflector.  Finally,  the  antenna  may 
cause  some  attenuation  due  to  the  gain  pattern  associated  with  an  orthogonal  polarization. 
A  reflected  signal  will  typiccdly  reverse  to  left  hand  circular  polarized  (LHCP)  while  the 
LOS  signal  will  be  RHCP. 

The  received  GPS  signal  in  the  presence  of  multipath  can  be  expressed  as 
rmpit)  =  y/^aom(t  —  To)c{t  -  To)cos{2Trfot  +  0o) 

M-l 

^  o<m(t  -  To  -  aiTc)c(t  -  To  —  oiiTc)cos{2Trfot  +  ${)  +  n(t)  (2) 
i=l 

where  P  is  the  power  of  the  signal,  M  represents  the  number  of  signals  received,  a,  is  the 
attenuation  coefficient  of  the  i*^  signal,  c{t)  is  the  spreading  code,  m{t)  is  the  navigation 
message,  fo  is  the  carrier  frequency  in  Hz,  To  is  the  propagation  delay  of  the  direct  path 
signal,  a^Tc  is  the  propagation  delay  of  the  signal  from  (relative  to)  the  direct  path 
signal,  0i  is  the  phase  of  the  signal,  and  n{t)  is  the  received  noise  signal.  The  noise 
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term  can  be  expressed  as  [9] 


n{t)  =  y/2  [nj  (t)  cos  (27r/ot)  -  ng  (i)  sin  (27r/ot)]  (3) 

where  ni{t)  and  nQ{t)  are  the  in-phase  and  quadrature  noise  components,  respectively. 

For  the  case  where  the  direct  path  signal  and  one  reflected  signal  is  received,  the 
signal  can  be  expressed  as 

rmp{t)  =  V^aom(t  -  To)c(t  -  To)cos(2Trfot  +  ^o) 

+'/^aim{t  —  To  —  aTc)c(t  — To  —  aTc)cos{2irfot  + 0i)  +  n{t)  (4) 

where  aTc  {a  =  ai  for  Equation  2)  represents  the  delay  of  the  reflected  signal  relative  to 
the  direct  path  signal.  The  average  power-delay  for  the  direct  path  signal  and  one  reflected 
signal  can  be  modeled  as  PaQ6{t  —  To)  -I-  Pai6{t  —  To  —  aTc)  as  shown  in  Figure  2. 


P 


Propagation  Delay 

Figure  2  Average  power-delay  for  direct  signal  and  one  reflection 

As  stated  earlier,  for  traditional  code  tracking  loops,  multipath  signals  with  delays 
less  than  l.STc  can  introduce  significant  errors  into  the  ranging  measurements.  For  ex¬ 
ample,  a  code  phase  tracking  error,  e,  of  one-tenth  of  a  chip,  O.lTc,  results  in  a  tracking 
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difference  of 


e  =  O.lTc  98ns 


(5) 


which  equates  to  a  ranging  error  of 

PRe  =  (3e8  m/s)(9.8e9  m/s) 

=  29.4  m.  (6) 


1.3.3  Receiver  Operation.  All  GPS  receivers  incorporate  some  sort  of  correlation 
process  in  their  design  in  order  to  remove  the  direct  sequence  spread  spectrum  (DS/SS) 
code  from  the  signal  (i.e., “despread”  the  signal).  This  receiver  despreading  operation  re¬ 
quires  a  correlation  of  the  received  signal  with  a  locally  generated  spreading  code  replica, 
which  uniquely  corresponds  to  the  desired  signal.  Conceptually,  a  received  signal  that  has 
been  spread  using  a  different  spreading  code,  corresponding  to  a  different  transmitter,  will 
appear  as  noise  and  will  cause  minimal  interference  with  the  desired  signal;  this  interfer¬ 
ence  corresponds  to  the  cross-correlation  between  the  two  spreading  codes,  which  can  be 
expressed  as 

rNTc 

Rcir)  =  1/NTc  /  c'  (t)  cit  +  T)  dt  (7) 

Jo 


where  N  represents  the  number  of  chips  per  code  cycle  and  both  waveforms  have  the  same 
chip  period  Tc. 

When  c(t)  =  (/(t),  a  sufficient  approximation  of  the  PRN  code  cross-correlation 
fimction  is  given  by  the  following: 


Rcir) 


^  1-^  for|r|<rc 
0  otherwise. 


(8) 


Observe  that  maximum  correlation  occurs  when  r  =  0,  corresponding  to  [perfect  synchro¬ 
nization.  Figure  3  illustrates  the  correlation  function  approximation  given  in  Equation  8. 
Note  that  the  peak  of  the  function  represents  maximum  correlation  between  the  locally 


6 


generated  spreading  code  and  the  received  signal.  As  the  two  signals  separate,  the  cor¬ 
relation  function  approaches  zero,  indicating  the  two  signals  are  no  longer  correlated.  As 
such,  the  received  signal  can  no  longer  be  detected  in  the  integration  of  Equation  7. 


Figure  3  Direct  sequence/ spread  spectrum  (DS/SS)  spreading  code  correlation  function. 


Previous  Research 

Currently,  errors  in  receivers  have  been  reduced  to  near  theoretical  limits  for  receivers 
processing  ideal  GPS  signals.  However,  non-ideal  signals  possess  characteristics  that  reduce 
positioning  accuracy.  These  errors  include  selective  availability  (SA)^,  satellite  clock  bias, 
atmospheric  effects,  ephemeris  errors,  and  multipath.  As  stated  previously,  the  influence 
of  most  of  these  errors  can  be  greatly  reduced  or  removed  entirely  through  differencing 
techniques.  Such  techniques  do  not  effectively  reduce  multipath  errors  because  of  the  lack 
of  multipath  correlation  between  receivers. 

Numerous  techniques  have  been  investigated  in  the  attempt  to  reduce  the  effects 
of  multipath.  Several  techniques  attempt  to  take  advantage  of  geometrical  differences 
between  the  direct  and  the  reflected  signals.  Other  techniques  involve  modifying  and 
creating  different  receiver  designs  by  enhancing  signal  processing  techniques.  This  thesis 
limits  its  analysis  to  techniques  associated  with  the  latter. 

The  delay-lock  loop  (DLL)  is  the  most  common  code  tracking  loop  design  used  in 
GPS  receivers.  It  is  composed  of  two  correlators  using  a  relative  delay  spacing  of  ATc 

^The  intentional  degradation  of  the  satellite  transmitted  signal  is  termed  selective  availability. 
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between  each  other.  The  ‘early’  correlator  mixes  a  code  replica  with  the  received  signal 
that  is  advanced  by  Arc/2  from  the  received  signal.  The  ’late’  correlator  mixes  a  code 
replica  with  the  received  signal  that  is  delayed  by  ATc/2  from  the  received  signal.  The 
signals  in  each  of  the  correlator  branches  are  processed  to  form  a  tracking  error  signal. 
Depending  on  how  the  correlated  signals  are  processed,  the  DLL  can  be  classified  as  being 
either  coherent  or  noncoherent.  The  coherent  DLL  requires  carrier  phase  information  to 
process  the  received  signals  whereas  the  noncoherent  DLL  (NCDLL)  operates  independent 
of  carrier  phase  tracking.  Traditionally,  receivers  employing  the  NCDLL  have  used  a 
normalized  correlator  spacing,  A  =  1.0.®  It  has  been  determined,  however,  that  narrowing 
the  correlator  spacing  improves  tracking  performance  [8,18]. 

Various  designs  have  been  derived  from  the  basic  properties  of  the  NCDLL.  However, 
new  designs  typically  include  estimation  hardware  to  determine  the  amplitude,  phase,  and 
delay  of  the  multipath  components  of  a  received  signal.  Two  such  designs  are  the  Multipath 
Estimating  DLL  (MBDLL)  and  the  Modified  RAKE  DLL  (MRDLL). 

MEDLL  estimates  the  amplitude,  delay,  and  phase  of  each  multipath  component 
using  ML  criteria.  The  estimates  are  used  to  remove  the  interfering  correlation  functions 
of  the  reflected  signals.  Thus,  only  the  direct-path  correlation  function  is  tracked  [14,16,17]. 

MRDLL  also  uses  ML  criteria  to  estimate  the  multipath  parameters.  These  estimates 
are  provided  to  a  multiple  correlator  tracking  loop.  Although  the  parameters  of  multiple 
reflected  signals  can  be  estimated,  MRDLL  has  only  been  simulated  imder  the  conditions 
with  one  reflected  signal  present  [6,7]. 

1.5  Scope 

This  research  investigates  the  use  of  ML  estimation  to  improve  code  tracking  per¬ 
formance.  MRDLL  serves  as  the  platform  in  which  the  ML  estimator  (MLE),  called  the 
multiple  correlator  estimation  unit  (MCEU),  is  implemented  and  tested.  In  addition,  the 
tracking  performance  of  MRDLL  is  compared  with  the  tracking  performance  of  the  NCDLL 
(with  a  normalized  correlator  spacing,  A  =  1.0  and  A  =  0.1). 

® ‘Normalized’  implies  a  spacing  A  relative  to  one  chip. 
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1.6  Assumptions 


The  following  assumptions  are  made  in  the  performance  of  this  research; 

1.  The  received  GPS  signal  consists  of  a  direct-path  signal  and  one  multipath  signal.^ 

2.  The  code  tracking  loop  has  acquired  the  phase  of  the  incoming  PN  sequence. 

3.  Only  C/A  code  is  being  tracked. 

4.  The  received  GPS  signal  corresponds  to  only  one  GPS  satellite. 

5.  Only  Signal  processing  techniques  are  considered  for  multipath  mitigation. 

6.  BPSK  data  modulation  on  the  received  signal  is  not  addressed. 

7.  Doppler  effects  on  code  treicking  performjmce  are  considered  negligible. 

1.7  Approach  Methodology 

An  in-depth  analysis  of  the  MCEU  is  performed  under  ideal  theoretical  conditions. 
An  enhanced  model  of  MRDLL  (referred  to  from  here  on  as  eMRDLL),  which  implements 
the  MCEU,  is  provided.  Using  the  model,  further  analysis  is  conducted  in  a  simulated 
code  tracking  environment  under  realistic  multipath  conditions.  In  addition,  the  track¬ 
ing  performance  of  eMRDLL  in  the  presence  of  additive  white  Gaussian  noise  (AWGN) 
environment  is  examined.  A  model  of  the  NCDLL  is  provided  to  investigate  the  track¬ 
ing  performance  of  the  NCDLL  (to  include  narrow  correlator  spa.cing)  versus  the  tracking 
performance  of  eMRDLL  in  a  simulated  multipath  environment. 

1.8  Materials  and  Equipment 

Theoretical  models  of  the  MCEU  were  implemented  using  Matlab,  version  4.2c,  and 
the  Optimization  Toolbox  from  The  MathWorks,  Inc.  of  Natick,  Massachusetts.  Sim¬ 
ulations  of  the  multipath  environment  and  models  of  eMRDLL  and  the  NCDLL  were 
developed  using  the  Simulink  Toolbox,  version  1.3,  the  Signal  Processing  Toolbox,  and 
the  Communications  Toolbox  which  are  also  from  The  MathWorks,  Inc.  Simulations  were 
performed  on  the  Sun  Workstations  provided  at  the  Air  Force  Institute  of  Technology 
(AFIT). 


■‘Note:  The  MLE  implemented  is  designed  to  estimate  the  parameters  of  multiple  reflected  signals. 


IL  Background 


2.1  Overview 

As  stated  in  Chapter  1,  numerous  techniques  have  been  applied  in  the  attempt  to 
reduce  the  effects  of  multipath.  Some  techniques  try  to  taJce  advantage  of  the  geometry 
of  the  direct  signal  and  the  reflected  signals,  by  varying  antenna  design  and/or  placement. 
Other  techniques  involve  modifying  and  creating  different  receiver  designs.  This  chapter 
examines  receiver  designs  that  incorporate  signal  processing  techniques  to  mitigate  the 
effects  of  multipath.  Improved  signal  processing  techniques  in  receiver  designs  offer  the 
greatest  potential  for  mitigating  the  effects  of  multipath  in  all  applications,  particularly 
with  the  recent  advances  in  technology. 

2.2  Delay  Lock  Loop 

One  of  the  predominantly  used  code  trax;king  loop  conflgurations  is  the  delay  lock  loop 
(DLL).  Commonly  referred  to  as  the  ‘early-late’  delay  lock  loop,  this  configuration  can  be 
designed  as  a  coherent  or  noncoherent  DLL  (NCDLL)  depending  on  the  desired  application. 
The  coherent  DLL  requires  a.ccurate  carrier  phase  tracking  whereas  the  NCDLL  can  operate 
without  carrier  aiding.  As  such,  the  NCDLL  is  less  likely  to  loose  lock  due  to  cycle  slips 
of  the  carrier,  thus  making  it  more  robust. 

Because  of  its  robustness,  the  NCDLL  is  employed  in  many  receiver  designs  to  per¬ 
form  code  loop  tracking.  It  is  composed  of  two  correlators  with  a  spacing  of  the  duration  of 
ATc  as  shown  in  Figure  4.  The  PN  generator  provides  the  ‘early’  gate  with  a  code  replica 
that  is  advanced  by  ATc/2  and  provides  the  ‘late’  gate  with  a  code  replica  that  is  delayed 
by  ATc/2  relative  to  the  received  signal.  Prior  to  narrow  correlator  spacing  advancements, 
a  correlator  spacing  of  Tq,  where  A  =  1.0,  was  typically  used.  The  correlator  outputs 
of  each  channel  are  bandpass  filtered  and  then  squared  and  low  pass  filtered  to  remove 
the  effects  of  data  modulation  and  carrier  phase  shift.  The  outputs  of  the  two  channels 
are  then  differenced  to  generate  the  code  phase  error  signal,  e{t,6),  which  is  subsequently 
filtered;  this  signal  is  used  to  drive  the  voltage  control  clock  (VCC)  which  corrects  the  code 
phase  error  of  the  local  PN  code  generator  (i.e.,  provides  feedback). 
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Figure  4  Non-coherent  Delay  Lock  Loop. 


2.2.1  Analysis  of  the  NCDLL.  The  received  signal  entering  the  NCDLL,  neglect¬ 
ing  the  navigation  message  and  multipath,  can  be  expressed  as 

r(t)  =  V^aoc{t  -  ro)cos{2irfot  +  ^)  +  n{t)  (9) 

where  n{t)  is  the  noise  term  presented  in  Equation  3.  For  this  analysis,  the  navigation 
message  can  be  neglected  because  this  results  in  the  maximum  possible  noise  component 
in  the  power  spectral  density  of  the  discriminator,  iSe,  within  the  tracking  loop  bandwidth. 
This  can  be  verified  by  realizing  that  without  data  modulation,  the  PSD  of  the  data, 
S0^{f),  is  a  delta  dirac  function  at  /  =  0  and  the  convolution  process  does  no  spectral 
widening  of  the  (signal  x  noise)  term  [9].  The  received  signal  enters  into  the  ‘early’  and 
‘late’  branches  of  the  NCDLL  where  it  is  mixed  with  a  locally  generated  replica  of  the 
spreading  code  represented  by  c(t  —  f©  ±  Arc/2)  in  which  the  +  sign  indicates  the  early 
and  —  sign  the  late  bremch  and  f©  represents  the  estimate  of  r©.  The  resulting  signal(s)  of 
the  ‘early’  and  ‘late’  breinches  are 

Xe{t)  =  V^Fooc(i  -  r©)c(t  —  f©  -h  ATc/2)cos{2irfot  +  (^)  +  n{t)  (10) 

and, 

xi{t)  =  V^aocit  —  To)c(t  -To  —  ATc/2)cos(2irfot  +  (f>)  +  n{t).  (11) 
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The  mixed  signals  in  each  channel  pass  through  a  bandpass  filter  (BPF)  centered 
at  /o  Hz.  The  BPF  cheuracteristics  should  be  chosen  such  that  it  passes  only  the  carrier 
frequency  plus  or  minus  any  Doppler  offset.  Simon  [11]  has  shown  that  code  tracking 
performance  may  be  improved,  particularly  in  deep  noise,  by  selecting  the  filter  bandwidth 
on  the  order  of  the  data  rate;  however,  Doppler  offset  must  still  be  taken  into  consideration. 

For  systems  with  high  spread  spectrum  processing  gain,  such  as  GPS,  the  components 
of  interest  are 

2/e(t)  =  V^aoc{t  -  To)c{t  -  To  +  ATcf2)cos{2Trfot  +  <!>)  +  ne{t)  (12) 


and 


yi(t)  =  v^Paoc(t  -  To)c(t  —  To  —  ATcf2)cos{2‘Kfot  +  0)  +  n/(t)  (13) 


where  the  noise  components  of  the  early  and  late  channels  are  given  by 


^e(^)  — 


c(t  -  To  +  -^Tc)n{t) 


(14) 


BPF 


and 


n/(t)  = 


c{t  -To  -  -;rTc)n{t) 


JBPF 


(15) 


respectively.  Although  code  self  noise  components  are  present  as  a  result  of  the  mixing 
operation,  these  components  can  be  neglected  for  systems  with  high  processing  gains  [9] 
and  where  the  single-sided  loop  bandwidth  is  much  less  than  the  PN  code  chip  rate 
(1.023  MHz  for  GPS  C/A  code)  [11]. 

Considering  that  the  dc  component  of  the  spreading  waveform  product  is  the  auto¬ 
correlation  function  Rdr)  of  the  spreading  waveform  evaluated  at  t  =  r©  —  To  ±  ATc/2, 

^GPS  loop  bandwidths  range  from  0.02  to  1  Hz  [19] 
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the  average  ‘early’  and  ‘late’  channel  signals  can  be  represented  as 


ye{t)  =  VPaoRc[{S  +  A/2)Tc]cos{2irfot  +  <!>)  +  ne{t)  (16) 

and 


Ut)  =  VPaoRc[{S  -  A/2)Tc]cos{2Trfot  +  (f>)  +  n7(t)  (17) 

where  6  =  (r  —  To)/3c  is  the  normalized  code  phase  estimation  error  and  n^{t),  ni{t)  denote 
time  averaged  noises.  The  signal  in  each  channel  passes  through  a  square  law  device  and 
then  is  low  pass  filtered.  The  LPF  removes  the  high  frequency  term,  47r/o,  caused  by 
the  squaring  operation  in  eax;h  channel.  The  difference  of  the  two  channels  provides  the 
tracking  discriminator/error  signal,  e(t,5)  =  [yfjiPF  —  [yt]LPF  given  as 

e(t,6)  =  Pa§5A(5)  +  ne(t).  (18) 


where  Sa(^)  is  the  delay  lock  loop  discriminator  (referred  to  as  the  S-curve)  defined  as 


Sa(S)  =  R^, 


(19) 


and  ne(t)  represents  the  noise  component  of  €(t, 5).  The  plot  of  the  tracking  error  for 
varying  values  t  forms  a  tracking  curve.  Different  values  of  ‘early-late’  gate  spacing  provide 
different  curve  characteristics,  such  as  a  chemge  in  slope  of  the  tracking  curve  which  in 
turn  causes  a  change  in  the  effective  tracking  range.  The  tracking  curve  is  linear  for  the 
normalized  correlator  spacing  A  <  1.0;  however,  it  becomes  non-linear  for  values  of  A  > 
1.0.  Therefore,  correlator  spacings  greater  than  1.0  chips  are  not  typically  implemented. 
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The  delay  lock  discriminator  for  A  <  1.0  can  be  characterized  as  [9] 

0  for— iV  +  H-^<5<  —(1  +  y) 

^  -  [1  +  (5  +  y)(1  +  ;^)]^  for  -  (1  +  y)  <  5  <  (y  -  1) 

-2{1  +  ^)A[1  + (1  +  ^)6]  for(f-l)<5<-f 

oaCo)  =  < 

2(1  +  i)[2  -  (1  +  ^)A]5  for  -  f  <  5  <  +1 

2(1  +  ^)A[1  -  (1  +  for  f  <  6  <  (1  -  f ) 

,  [1  -  (1  +  -  t)]^  -W  for  (1  -  f )  <  ^  <  (1  +  t) 


For  the  traditional  value  of  A  =  1.0,  the  resulting  tracking  curve  is  presented  in  Figure  5. 


Figure  5  S-curve  for  NCDLL,  A  =  1 


2.2.2  Determination  of  Noise,  n^it).  The  noise  components,  ne(t)  and  ni{t),  at 
the  output  of  the  BPF  are  the  products  of  a  bandlimited  white  Gaussian  noise  process  n(t) 
centered  at  /©  and  the  spreading  code  waveforms  c(t— foiyTc)  as  described  in  Equations  14 
and  15.  This  process  provides  the  benefit  of  spreading  the  noise  via  the  spreading  code 
waveforms  and  thereby  reducing  the  power  spectral  density  in  the  frequency  range  passed 
by  the  BPF.  The  power  spectrum  of  n{t)  after  passing  through  the  BPF  ceui  be  expressed 
as 

Sn':eAf)  =  Sn:eAf)mf)\^-  (21) 
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Since  the  power  spectrum  of  n(t)  is  much  wider  in  bandwidth  than  the  power  spectrum 
of  c{t),  convolution  has  little  effect.  Assuming  ideal  filters,  the  output  PSD  of  the  noise 
component,  Sn>:e,l{f),  has  a  magnitude  of  ^  and  two  sided  bandwidth  Bn  and  is  centered 
about  /o. 

The  noise  at  the  discriminator  output,  ne(t),  is  a  function  of  e(t,6)ng  =  [yf]LPF  — 
[yg]LPF'  As  stated  e2irlier,  when  the  received  carrier  is  not  modulated  with  data,  the 
maximum  possible  noise  component  is  produced  in  the  power  spectral  density,  5^,/), 
within  the  tracking  loop  bandwidth  [9].  Because  the  loop  filter  bandwidth,  Bi,  following 
the  discriminator  is  much  smaller  than  ne(t),  the  dc  component  of  5^,(7)  provides  an 
accurate  approximation  of  the  noise  PSD.  This  PSD  is  approximated  as  [4,9] 


SnAf)«Snmf=0  =  2NiBN 

+2PagjVo{iJ^[(«-|)re  -■Ro[(^  +  f)rc]}  (22) 


over  the  loop  bandwidth  (zero  elsewhere)  and  has  units  of  W/Hz. 


2.2.3  NCDLL  Linear  Model  The  delay  lock  discriminator  is  determined  by  a 
non-lineax  process  (because  of  the  squaring  operation  in  the  square  law  device).  However, 
a  linear  model,  which  will  be  developed  in  this  section,  can  still  be  used  to  describe  the 
DLL. 

Ultimately,  the  code  phase  tracking  error  is  determined  by  the  difference  of  the  phase 
of  the  received  signal  and  the  code  phase,  produced  by  the  VCC,  and  is  given  as 


Tdit)  =  Ki6 

=  if.(^)  P8) 

where  S  represents  the  normalized  code  phase  error,  To/Tc  represents  the  normalized  code 
phase  of  the  received  signal,  fo /Tc  represents  the  normalized  estimated  code  phase  from 
the  VCC,  and  Kd  represents  the  gain  of  the  mixing  process  in  volts.  The  VCC  operates 
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at  an  angular  frequency  represented  as 


27r/vc'C'(*)  =  27r/,  +  Korf{t) 


(24) 


where  2'Kfq  is  the  center  frequency  (rad)  of  the  VCC,  Ko  is  the  VCC  gain  in  {rad/sec  volts), 
and  r/(t)  is  the  output  signal  of  the  loop  filter  [1].  When  6  =  0,  the  VCO  is  operating  at 
the  center  or  quiescent  frequency,  27r/g.  If  the  input  signal  is  positive  at  a  specific  instant 
in  time,  the  phase  of  the  input  signal  then  starts  leading  the  phase  of  the  output  signal  of 
the  VCC.  An  error  signal,  e(t,  6),  develops  at  the  discriminator  output  which  increases  with 
time.  With  a  delay  due  to  the  loop  filter,  r/(i)  will  increase  with  time  causing  the  VCC 
to  increase  its  frequency  until  the  VCC’s  frequency  matches  the  input’s  signal’s  frequency. 
The  VCC  then  operates  at  a  new  frequency  equal  to  27r(/o + A/).  The  signed  at  the  output 
of  the  loop  filter,  r/(t),  will  settle  to  a  final  value  of  rf  =  2irAf/Ko. 

The  angular  frequency  of  a  signal  is  defined  as  the  first  derivative  of  its  phase  with 
respect  to  time,  2Trf  =  d$fdt.  Hence,  the  phase  of  a  signal  is  the  time  integral  of  its 
angular  frequency.  Equation  24  integrated  over  time  is  [9] 


which  simplifies  to 


To 


27r;jfr(<)  = 


I  Tf{pi)doc 

Jo 


To(t) 


=  Ko  I  rf{a)da. 

Jo 


(25) 


(26) 


where  fo{t)/Tc  represents  the  phase  of  the  VCC.  The  output  of  the  loop  filter  can  be 
described  by  its  impulse  response  f{t)  convolved  with  the  input  signal  as  follows 


r/(a)  =  e{t,S)  *  f{X)d\ 


Incorporating  Equation  27  into  Equation  26,  the  output  phase  of  the  VCC  cam  be  expressed 
as 


=  Ko  f  f  e{a,6)f{a— \)dXda. 
Jo  J—oo 


Substituting  Equation  18  into  Equation  28  provides  the  equation  for  the  nonlinear  model 


Ko  f  [  {PalSA{S)  +  ne(A))  /(a  -  X)d\da. 

Jo  J — oo 


For  small  code  tracking  errors,  —4  <  5  <  +4,  Equation  20  provides 


S.W  =  2(i  +  1)[2-(i+1)a]« 


To  form  a  linear  model,  shown  in  Figure  6,  the  noise  term,  ne(t),  of  the  discriminator 
is  moved  to  the  loop  input  with  the  appropriate  gain  adjustment  provided  as  follows 

=  +  +  (31) 

Based  on  the  linear  model  shown  in  Figure  6,  the  Laplace  transform  of  the  tracking  loop 


Figure  6  NCDLL  linear  equivalent  circuit. 


output  To  relative  to  the  loop’s  input,  To,  can  be  represented  via 


To(s)  ro(s)  -  fo(s) 


KdK, 


F(s) 


(32) 


Solving  for  given  Equation  32  results  in  the  closed  loop  transfer  function 

.  A  fo(s)  KdKoFjs) 

To{s)  8  +  KdKoF{s) 


(33) 


which  characterizes  the  DLL.  Commonly,  a  simple  first  order  lead-lag  filter  is  implemented 
as  the  loop-filter,  F{s).  The  transfer  function  of  the  loop  filter  can  be  expressed  as 


Fis)  = 


1  +  T2S 
ai  +  Tis 


(34) 


where  [3] 


ai 


■{: 


passive  filter 
0  active  filter 


(35) 


It  has  been  determined  that  better  performance  will  almost  always  be  obtained  through 
the  use  of  an  active  filter  [2].  Assuming  the  use  of  a  first-order  active  filter,  produces  a 
second  order  loop  transfer  function 


= 


(^) 

KaKo  ,  K^KoV, 
ns  n 


ns 


n 


I  K^Ko 


(36) 


Equation  36  can  then  be  represented  in  the  classical  manner  for  servo  mechanisms  as  [2,9] 

2Cu;nS  + 


His)  = 


s2  +  2Cfa;wS  + 


(37) 
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where  the  loop  natural  frequency,  in  and  the  damping  factor,  ^  (unitless),  are 
defined  as 


and 


c  = 


(39) 


respectively.  The  single-sided  noise  equivalent  bandwidth,  Bl  in  Hz,  for  a  second-order 
loop  with  an  active  lead-lag  loop  filter  can  be  described  as  [2,8] 


Bl 


poo 

/  \Bij27rf)fdf 
Jo 


(40) 


The  natural  firequency,  u;„,  and  the  damping  factor,  determine  the  tracking  performance 
of  the  NCDLL.  Increasing  the  natural  frequency  will  decrecise  the  loop  response  time.  The 
noise  loop  bandwidth,  Bl,  determines  the  amount  of  noise  power  that  effects  the  tracking 
loop.  Therefore,  loop  bandwidth,  Bl,  is  chosen  by  the  designer  so  as  to  minimize  the 
effects  of  noise  balanced  with  adequate  response  time.  Smaller  loop  bandwidths  reduce 
noise  levels,  while  increasing  response  time  resulting  in  more  sluggish  responses. 


2.2.4  NCDLL  Tracking  Performance  Characteristics.  Mean-square  tracking  er¬ 
ror  or  tracking  jitter  is  of  particular  interest  in  ainalyzing  the  performance  of  the  NCDLL. 
The  power  spectrum  of  the  variance  of  the  steady  state  tracking  error,  S,  is  given  by  [2,9] 


SsU)  =  l^(i27r/)|25„,-(/) 


(41) 


where  H (s) \t=j2irf  is  the  closed-loop  transfer  function  defined  in  Equation  37  and  Sne'if)  = 
5n«(/)/iirJ  is  the  two-sided  power  spectrum  of  the  Gaussian  noise  process  at  the  input  to 
the  loop  model.  Ss{f)  is  approximately  flat  over  the  loop  bandwidth,  Bl,  so  that  the 
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variance  of  8  can  be  expressed  as 


/oo 

Ss{f)df 

‘OO 

=  r  SnAf)mj2irf)\‘^df 

«/— OO 

where  the  integral  on  the  right  side  of  the  equation  is  defined  as  the  two-sided  noise 
bandwidth  Bi  in  Hz.  For  A  <  1.0,  the  variance,  known  as  tracking  jitter,  then  has  the 
form  [8,18] 


BlA 

2SINq 


1  + 


2 

(2  -  A)5/(J\roBjv). 


(43) 


where  Bu  represents  the  single-sided  bandwidth  of  the  BPFs  and  LPFs.  Equation  43 
indicates  that  eis  A  decreases,  the  tracking  jitter  variance  decreases.  This  is  caused  by  the 
cancelation  of  noise  due  to  the  overlap  in  the  early  and  late  channels  for  A  <  1.0. 


2.2.5  Performance  of  the  NCDLL  With  Multipath.  The  NCDLL  performs  well 
in  tracking  the  direct  path  signal  with  no  multipath.  However,  when  multipath  delays 
are  present  within  a  range  of  (0,1.5]  chips,  tracking  errors  may  become  significant.  This 
section  examines  the  effects  of  multipath  on  NCDLL  performance.  It  is  assumed  that  only 
one  reflected  signal  is  present.  Neglecting  the  navigation  message,  as  done  previously,  the 
received  signal  with  one  reflection  can  be  expressed  eis 

'I'mpi't)  =  V^aocft  —  To)cos(27r/o<  +  ^o) 

-|-V^Poic(t  —  To  —  aTc)cos{2Trfot  9i)  -f-  n{t).  (44) 


where  the  second  term  represents  the  reflected  signal  with  delay  aTc.  As  in  the  previous 
section,  the  received  signal  enters  into  the  ‘early’  and  ‘late’  branches  of  the  NCDLL  upon 
which  it  is  mixed  with  a  locally  generated  replica  of  the  spreading  code  represented  by 
c{t  —  To  i  Arc/2)  in  which  the  -I-  sign  indicates  the  ‘early’  branch  and  —  sign  the  ‘late’ 
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branch  giving  the  resulting  signals 


Xe(t)  =  \/^aoc(t  —  ro)c(t  —  To  +  ATc/2)cos(27rfot  + 

+V^aic(t  —  To  —  aTc)c(t  -To  +  A.Tc/2)cos{2Trfot  +  (f>)  +  ne(t)  (45) 


and, 

xi(t)  =  y/^aoc(t  —  To)c(t  —  To—  ATc/2)cos(2Tfot  +  (f>) 

+y/Waic{t  —  To  —  ocTc)c{t  -To  —  ATc/2)cos{2nfot  +  (j>)  +  n/(t).  (46) 


Once  again,  for  high  DS/SS  processing  gain  systems  such  as  GPS,  code  self  noise  can 
safely  be  neglected.  Therefore,  the  signals  after  bandpass  filtering  becomes 

ye{t)  =  V^aoRc  ^cj  cos{2Trfot  +  6o) 

+VWaiRc  +  a  +  Tcj  cos{2nfot  +  Bi)  +  ne(t)  (47) 

and, 

yi{t)  =  V^aoRc  Tcj  cos{2nfot  +  ^o) 

+V^aiRc  +  a  -  y^  Tcj  cos(27r/ot  +  6i)  +  ni(t)  (48) 

where  S  =  {t  —  To)ITc.  The  squaring  and  LPF  operation  produces  three  signal  components 
for  each  branch  which  is  a  result  of  squaring  the  direct  path  signal  and  the  reflected  signal 
and  a  cross  term  of  the  direct  path  and  reflected  signal.  Neglecting  the  noise  terms,  the 
signals  can  be  expressed  as 


yl 


p4ri]^I^6+j^Tc^ 

‘\-2PaQaicos{0o  —  0i)Rc  [(«  +  f)r.]je.[(«  +  c,+  f)r,] 

+a5i^  [(«  +  «+ |)n] 


(49) 
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and 


2/?  = 


PalRl  («-f)rc] 

+2Paoaico8{6o  —  6i)Rc  Tcj  Rc  +  a  —  Tcj 

+a2ij2  [(«  +  «- y)r,] 


(50) 


for  the  ‘early’  and  ‘late’  branches  respectively.  Differencing  the  two  branches  forms  the  S 
curve,  €(<,  S)  =  yf  —  which  can  be  expressed  as 


[(«  +  !)  r,]  I 

-h2Pa[iaia>s(0o  —  ^l) 

+  [(«  +  «  + f)r.]} 

+Paf  [(s  +  a  -  ^)  r,]  -  [(i  +  a  +  I)  T.]  j .  (61) 

where  ffo  —  0i  =  2irfoaTc.  The  first  and  third  terms  represent  the  direct-path  and  multipath 
components,  respectively.  The  second  term  is  the  cross  term  component  caused  by  the 
interaction  of  the  direct-path  and  multipath  signal.  Figure  7  illustrates  the  influence  of 
the  multipath.  It  should  be  noted  that  the  direct-path  component  is  centered  about  5  =  0. 
However,  the  net  influence  of  the  other  components  generated  as  a  result  of  multipath, 
generates  a  tracking  error,  6*  ^  0,  corresponding  to  e(5*)  =  0.^  The  net  result  of  the 
three  components  is  illustrated  in  Figure  8.  Here,  one  can  see  that  the  tracking  point, 
corresponding  to  €(5*)  =  0,  is  no  longer  at  5  =  0,  but  has  some  multipath  induced  tracking 
error.  The  result  is  a  steady-state  code  pheise  tracking  error,  6*,  where  6*  =  .  The 

amplitude  of  the  code  phase  tracking  error  depends  on  the  amplitude  and  phase  parameters 
of  the  direct-path  and  multipath  components:  uq,  ai,  0o>  and 

The  predicted  steady-state  tracking  error  can  be  determined  by  setting  Equation  51 
equal  to  zero  and  solving  for  6  for  any  a  6  [0,1.5].  Figure  9  provides  the  predicted  steady- 

^6*  lepresents  the  actual  tracking  point  as  opposed  to  the  desired  tracking  point. 

® Recall  that  a2itft,Tc  =  6o  —0i. 


■(6-f)r.]  A,  [(«  +  »- f)?; 
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Figure  7  Typical  NCDLL  discriminator  output  components  in  the  presence  of  multipath. 

state  tracking  error^  for  ao  =  1.0,  ai  =  0.5,  P  =  0.5,  and  foTc  —  10.  As  presented  by  van 
Nee  (for  A  =  1.0),  multipath  delays  of  up  to  1.5  chips  can  cause  errors  in  the  tracking 
loop  estimates  [19].  The  figure  also  illustrates  that  maximum  tracking  errors  occur  when 
a  multipath  signal  is  in-phase  or  out-of-phase  with  the  line  of  sight  signal,  “  ^1  =  Trn, 
where  n  is  any  integer.  The  estimate  of  the  carrier  phase,  however,  will  theoretically  have 
no  tracking  error,  because  carrier  tracking  follows  the  combination  of  the  direct-path  and 
multipath  signal  which  have  an  equivalent  zero  crossing  point  when  they  are  in-phase  or 
180°  out  of  phase  [15,19].  GPS  has  a  product  of  foTc  =  1540.  Therefore,  the  number  of 
cycles  of  the  tracking  error  curve  would  be  much  higher,  within  the  same  tracking  error 
envelope;  the  envelope  around  the  tracking  error  curve  indicates  the  worst-case  (i.e., in- 
phase  or  out-of-phase  multipath  signal)  for  any  product  of  /oTc. 

2*2,6  Narrow  Correlator  Spacing,  One  of  the  early  advances  in  receiver  design 
to  combat  multipath  was  proposed  by  A.J.  Van  Dierendonck,  Pat  Fenton,  and  Tom  Ford 
in  the  early  1990’s  [18].  They  proposed  narrowing  the  chip  spacing  in  the  NCDLL.  This 
came  about  after  evaluating  the  traditional  design  of  the  NCDLL. 

^Note:  The  envelope  around  tracking  error  curve  indicates  worst-case  values  for  any  product  of  foTc, 


23 


Figure  8  Typical  NCDLL  discriminator  output  in  the  presence  of  multipath 

The  improvements  obtained  by  narrow  correlator  spacing  of  O.lTc  require  the  use 
of  a  wide-bandwidth  receiver  design  to  sharpen  the  peak  of  the  code  correlation  function. 
Due  to  the  close  spacing  of  the  reference  codes,  the  noise  in  the  ‘early’  and  ‘late’  correlator 
outputs  is  highly  correlated.  As  such,  the  noise  is  reduced  when  the  ‘early’  and  ‘late’  paths 
are  differenced.  In  addition,  in  the  presence  of  multipath,  the  error  tracking  envelope  is 
smaller  which  leads  to  improved  code  tracking  performance  [20]. 

In  the  traditional  design,  it  was  assumed  that  a  one  chip  spacing  provided  optimum 
performance.  In  the  previous  designs,  analog  hardware  was  used  in  GPS  receivers.  A  one 
chip  spacing  minimized  hardware  requirements.  Also,  early  receivers  were  designed  for 
receiving  P-code  which  has  a  short  chip  width.  Implementing  narrower  spacings  makes 
the  delay  lock  loop  discriminator  very  narrow.  It  was  feared  that  Doppler  shift  and  other 
disturbances  would  cause  a  loss  of  code  lock.  Neurower  spacing  eilso  requires  faster  clock 
rates  to  generate  code  replicas  for  the  early  and  late  gates.  Previously,  technology  limited 
the  ability  to  adjust  for  these  considerations.  Today’s  technology  allows  the  use  of  narrower 
chip  spacing  between  the  early  cuid  late  gate  for  C/A  code  [18].  As  shown  in  Figure  10, 
the  slope  of  the  tracking  curve  is  much  steeper  for  A  =  0.1  than  the  slope  for  A  =  1.0 
which  indicates  the  requirement  of  a  more  dynamic  VCC. 

A  closed  form  solution  of  the  tracking  envelope  for  the  coherent  DLL  has  been  derived 
which  can  be  used  to  examine  the  benefits  of  narrow  correlator  spacing.  The  solution 
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Figure  9  Predicted  NCDLL  multipath  tracking  error  for  parameters:  ao  =  1.0,  ai  =  0.5, 
P  =  0.5,  foTc  =  10,  and  A  =  1.0. 

can  be  used  to  analyze  the  NCDLL  as  well  since  the  worst-case  tracking  errors  for  a 
noncoherent  DLL  are  the  same  as  for  a  coherent  DLL  [15].  By  analyzing  the  correlation 
function  for  the  direct  path  signal  in  conjunction  with  the  reflected  signal,  four  sets  of 
equations  can  be  obtained,  summarized  in  Table  1,  through  algebraic  manipulation  [15]. 
These  equations  are  based  on  the  interaction  of  the  direct  path  correlation  function  with 
the  reflected  signals  correlation  function  through  four  regions.  The  equations  are  valid  for 

Table  1  Equations  for  determining  the  code  tracking  phase  error  envelope  of  a  coherent 
or  NCDLL. 


Range  of  a 

6 

JhaSL. 

an4~am 

i22+^<a<re-A(l-i2^)  (=6) 

am  A 

2ao 

a  >  Tc  H-  f 

0 

the  condition  that  A  <  Tc  euid  — ao  <  <  oq.  The  constant  corresponds  to  ±ai 

depending  upon  whether  the  reflected  signal  is  in  phase  or  out  of  phase  with  the  direct 
path  signal.  Figure  11  displays  the  points  of  the  tracking  error  envelope  generated  from 
Table  1  which,  in  conjimction  with  Figure  11,  can  be  used  to  describe  the  tracking  error 
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S-curve 


Figure  10  S-curve  for  NCDLL,  A  =  0.1 


envelope  for  any  single  multipath  reflection  with  any  veilue  of  A.  Figure  11  illustrates  the 
reduction  in  tracking  error  afforded  by  narrowing  the  correlator  spacing,  A. 


Figure  11 


Code  tracking  error  envelope  (worst-case)  in  terms  of  multipath  parameters 
and  A.  Note:  ±a  and  ±b  axe  the  points  relative  to  ±oi  for  the  equations  in 


Table  1. 


Figure  12  provides  the  results  of  narrowing  A  to  0.1  with  the  same  signal  parameters 
as  used  in  Figure  9.  Upon  comparison,  the  worst  case  tracking  error  is  reduced  by  an 
order  of  magnitude  and  the  range  of  a  which  produces  tracking  error  has  been  reduced  to 
(0,1.05].  Note  that  the  frequency  of  the  tracking  error  curve  is  the  same  as  in  Figure  9, 
foTc  =  10.  Figure  13  provides  an  overlay  of  the  tracking  error  envelopes  for  A  =  1.0 
and  A  =  0.1  which  clearly  indicates  the  significant  improvement  afforded  by  naxrowing 


the  correlator  spacing.  As  such,  the  narrow  correlator  design  of  the  NCDLL  provides  a 


Figure  12  Predicted  NCDLL  multipath  tracking  error  for  parameters:  cq  =  1.0,  ai  = 
0.5,  P  =  0.5,  foTc  =  10,  and  A  =  0.1. 


Figure  13  Code  tracking  error  envelopes  for  A  =  1.0  and  A  =  0.1 

very  good  baseline  for  evaluation  of  other  multipath  mitigation  techniques.  These  NCDLL 
properties  have  been  documented  by  many  sources  including  [8,9,12].  The  NCDLL  with 
narrow  correlator  spacing  will  be  used  as  a  performance  baseline  this  thesis  as  well. 

2.3  Estimator  Designs 

Although  the  NCDLL  with  narrow  correlator  spacing  has  had  a  profoimd  effect 
on  multipath  reduction,  better  receiver  designs  have  been  achieved  using  estimators  to 
determine  the  parameters  of  multipath  signals  in  an  effort  to  totally  remove  the  error 
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caused  by  multipath.  Three  designs  using  estimation  theory  to  remove  the  effects  of 
multipath  in  a  DS/SS  or  Code  Division  Multiple  Access  (CDMA)  environment  are  the 
multipath  estimating  DLL  (MEDLL),  RAKE  DLL  (RDLL),  and  the  Modified  RAKE  DLL 
(MRDLL). 

2.3.1  MEDLL.  The  MEDLL  was  introduced  by  Richard  van  Nee  to  improve  the 
performance  of  positioning  solutions  in  a  multipath  environment  [16,17].  This  design  is 
based  on  Maximum  Likelihood  Estimation  (MLE)  theory  which  is  used  to  estimate  the 
phase,  propagation  delay,  and  amplitude  of  the  direct  path  and  refiected  signals.  One 
method  of  implementing  MEDLL  is  accomplished  by  using  a  bank  of  M+1  correlators, 
where  M  is  the  number  of  reflected  signals,  to  simultaneously  estimate  the  parameters  of 
the  line-of-sight  signal  as  well  as  the  reflected  signals  [16].  The  MEDLL  determines  a  set 
of  reference  correlation  functions  with  a  certain  amplitude,  phase,  and  delay  which  gives 
the  best  estimate  of  the  input  correlation  function.  Eax;h  estimated  multipath  correlation 
function  is  subtracted  from  the  measured  correlation  function  which  provides  a  remaining 
estimate  of  the  direct  path  correlation  function.  A  standard  coherent  early-late  DLL 
proceeds  this  process  to  provide  an  estimate  of  the  code  loop  tracking  error.  Essentially, 
all  equations  look  the  same  as  in  the  case  of  a  conventional  coherent  spread  spectrum 
receiver  except  for  the  estimation  and  removal  of  multipath  signals  [14,16].  In  normal 
operation  the  MEDLL  is  configured  to  estimated  the  parameters  of  the  direct  path  signal 
and  two  reflected  signals  [14]. 

2.3.2  RAKE  Delay  Lock  Loop.  The  basic  configuration  of  the  DLL  was  not  in¬ 
tended  to  be  used  under  adverse  conditions  such  as  envelope  fading,  code  delay  spread,  and 
code  Doppler  spread.  Stuber  and  Sheen  [10, 13]  proposed  a  new  low  complexity  tracking 
loop  for  direct-sequence  spread-spectrum  signaling  on  multipath  fading  inter-symbol  inter¬ 
ference  (ISI)  channels  found  in  applications  such  as  code-division-multiple-eiccess  (CDMA) 
mobile  radio  communications.  The  new  tracking  loop,  called  the  RDLL,  exploits  the  inher¬ 
ent  multipath  diversity  of  the  channel  similar  to  a  RAKE  receiver.  Although  not  intended 
for  GPS,  the  RDLL  has  been  shown  to  have  excellent  performance  as  compared  with  the 
traditional  DLL  for  ranging  applications  such  as  mobile  radio  communications  [10,13]. 


28 


The  RDLL  consists  of  a  channel  parameter  estimation  unit  and  a  coherent  tracking 
loop.  The  channel  pareimeter  estimation  unit  is  composed  of  L+1  branches  spaced  at 
integer  multiples  of  Tc  which  provide  estimates  to  the  correlator  branches  after  mixing.  The 
channel  parameter  estimation  is  a  moving  average  filter  that  averages  the  received  signal 
following  despreading  and  demodulating.  The  best  averaging  length  can  be  empirically 
determined  and  is  based  on  the  Doppler  shift  range  of  values  of  the  channel  [13].  The 
coherent  tracking  loop  consists  of  L+1  correlator  functions  which  use  an  ‘early-late’  code 
correlation  process.  The  correlated  signals  are  converted  to  basebeind  euid  summed  to  form 
a  low  pass  error  signal.  This  signal  drives  a  Voltage  Control  Clock  (VCC)  which  corrects 
the  code  phase  error  of  the  locally  generated  pseudo-noise  code  [10, 13). 

2.3.3  Modified  RAKE  Delay  Lock  Loop.  The  MRDLL,  a  design  proposed  first 
proposed  by  Laxton  and  DeVilbiss,  uses  maximum  likelihood  signal  estimation  to  deter¬ 
mine  the  signal  parameters  of  the  direct  path  signal  and  a  single  reflected  signal  [6,7].  The 
MRDLL,  shown  in  Figiire  14  is  a  modified  version  of  the  RDLL  introduced  by  Sheen  and 
Stuber.  It  is  composed  of  three  main  components:  the  multiple-correlator  tracking  loop 
(MCTL),  the  adaptive  loop  controller  (ALC),  and  the  multiple  or  multipath  correlator 
estimation  unit  (MCEU). 


MCTL 


Figure  14  Modified  RAKE  delay  lock  loop  (MRDLL)  [6,7]. 
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After  conversion  to  baseband,  the  received  direct  path  signal  is  tracked  by  the  MCTL. 
The  MCTL  relies  on  the  MCEU  generated  maximum  likelihood  estimated  signal  param¬ 
eters  to  accurately  track  the  direct  path  code  phase  in  the  presence  of  multipath.  It  was 
conceived  as  a  bank  of  M  tracking  arms  of  which  one  is  designated  for  the  direct  path 
signal.  The  remaining  tracking  arms  are  delayed  by  some  range  of  values  from  the  direct 
path  trax:king  arm.  In  implementation  a  movable  eirm  was  actually  used.  The  discriminator 
output  is  formed  by  summing  the  direct  path  arm  with  the  arm  that  has  been  determined 
to  contain  the  most  accurate  estimate  of  the  correlation  function  of  the  reflected  signal. 
The  MCEU  is  composed  of  two  components:  a  bank  of  M  correlators  and  an  estimator. 
It  is  designed  to  provide  the  MCTL  with  maximum  likelihood  estimates  of  the  multipath 
delay  with  respect  to  the  direct  path  delay  as  well  as  coefiicients  of  the  LOS  baseband 
signal  and  the  reflected  baseband  signal.  These  estimates  allow  the  MCTL  to  remove  the 
tracking  error  introduced  by  the  reflected  signal.  Whereas  the  RDLL  models  the  relative 
delays  between  the  direct  path  and  reflected  signals  as  integer  multiplies  of  the  spreading 
code,  it  has  been  shown  that  smaller  delays  of  fractions  of  a  chip  cause  tracking  errors 
in  GPS  applications.  Therefore,  the  MCEU  incorporates  chip  spacings  in  a  range  of  (0, 
1.5].  The  MCEU  feeds  the  estimate  of  the  direct  path  signal  coefficient,  ®o>  to  the  direct 
path  arm  in  the  MCTL.  It  also  feeds  the  estimate  of  the  reflected  signed  coefficient,  ®i,  to 
the  movable  tracking  arm  which  is  then  used  to  provide  the  most  accurate  estimate  of  the 
correlation  function  of  the  reflected  signal.  Estimates  of  both  signal  coefficients  as  well  as 
the  estimated  delay  of  the  reflected  signal  Eire  provided  to  the  ALC. 

Enhancements  to  the  MRDLL  have  improved  its  performance.  The  enhanced  MRDLL 
(eMRDLL)  serves  as  the  basic  design  for  which  performance  of  maximum  likelihood  es¬ 
timation  techniques  can  be  analyzed.  The  methodology  for  the  eMRDLL  is  presented  in 
Chapter  3. 
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III.  Enhanced  Modified  RAKE  Delay  Lock  Loop  (eMRDLL) 


3.1  Overview 

The  enhanced  modified  RAKE  delay  lock  loop  (eMRDLL),  shown  in  Figure  15,  is 
composed  of  three  main  components:  the  multiple-correlator  tracking  loop  (MCTL),  the 
adaptive  loop  controller  (ALC),  and  the  multiple  or  multipath  correlator  estimation  unit 
(MCEU).  This  coherent  design  requires  the  received  signal  to  be  converted  to  baseband 
prior  to  or  during  the  correlation  processes  within  the  MCTL  and  the  MCEU. 

The  ‘early’  ‘late’  gate  tracking  loops  of  the  MCTL  employ  the  same  basic  charac¬ 
teristics  as  a  coherent  DLL  with  the  added  capability  of  correcting  the  tracking  process 
according  to  the  estimated  amplitudes  and  delays  of  the  received  LOS  and  multipath  sig¬ 
nals  via  the  MCEU.  The  ALC  dynamically  adjusts  the  loop  tracking  response  according 
to  the  estimated  signail  parameters. 
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3.2  Baseband  Conversion  Process 

The  first  processing  stage  of  eMRDLL  converts  the  received  multipath  signal  to 
baseband.  The  received  multipath  signal,  rmp{t),  can  be  expressed  as  follows 

Tmp{t)  =  V^aoc{t  -  To)cOs{2'Kfot  +  0o) 

+V^aic(t  —  To  —  aTc)cos(2Trfot  +  0i)  +  n{t)  (52) 


where  n(t)  represents  the  noise  expressed  in  Equation  3.  As  part  of  the  baseband  conversion 
process,  the  signal  is  mixed  with  a  locally  generated  signal,  2cos(27r/o  +  ^3),  produced  by  a 
phase  lock  loop  (PLL),  and  then  low  pass  filtered  in  the  MCTL  and  MCEU  to  convert  the 
multipath  signal  to  baseband.  This  represents  a  modification  in  the  baseband  conversion 
process  as  presented  by  Laxton  and  DeVilbiss,  which  contained  a  pre-correlation  filter  to 
convert  the  received  signal  to  baseband  prior  to  entering  the  MCTL  and  MCEU  as  shown  in 
Figure  14  in  Chapter  2  [6,7].  Using  a  pre-correlation  filter  is  not  necessary  since  the  mixed 
signal  will  be  low  pass  filtered  within  the  MCTL  and  MCEU.  In  addition,  low  pass  filtering 
prior  to  the  correlation  process  may  distort  the  signal,  r{t).  Therefore,  implementation  of 
a  pre-correlation  filter  would  require  the  coefficients  of  the  filter  to  be  considered  in  the 
estimation  process  in  order  for  the  estimator  in  the  MCEU  to  meet  the  Cramer  Rao  lower 
bound  (CRLB)  as  a  minimum  variance  unbiased  (MVU)  estimator  [22]. 


The  PLL  signal  can  be  generated  by  passing  the  input  signal,  rmp{t),  through  a 
correlator  followed  by  a  phase-locked  loop  (PLL)  as  shown  in  Figure  16  to  produce  the 
signal,  2cos(27r/oi  +  ^3)  [6]. 


Tmpij') 


BPF 


im 


PLL 


2cos{2'Kfot  +  9^) 

(to  baseband  mixer) 


c{t  -  To) 


Figure  16  Typical  carrier  phase  recovery  scheme  for  GPS. 
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The  signeil,  entering  the  PLL,  obtained  through  the  correlation  process  of  mix¬ 
ing  and  bandpass  filtering  the  received  signal,  can  be  expressed  as 

r'{t)  =  ^V^Faoc(t  -  ro)c(t  —  fo)cos(27r/ot  -f  ^o) 

-h  V^aic{t  -To  —  aTc)c{t  —  fo)co8{2irfot  +  ^i)] 

J  BPF 

=  y/^aoRc{0)cos{2Trfot  +  0o)  +  y/^aiRc{a)cos{2Trfot  +  6i) 

=  V^aocos(2Tfot  + 0o)  +  V^aiRc{oi)cos{2Tfot  +  $i)  (53) 

Assuming  perfect  code  phase  synchronization,  To  =  To,  an  expression  for  ^3  CcUi  be  derived 
by  the  following  technique  [23].  By  replacing  27r/oi  -b  0i  with  27r/ot  -b  ^0  +  ~  ^0  and  using 

the  trigonometric  identity  cos{u  ±v)  =  cos{u)cos{v)  ^  8in(u)sin{v),  Equation  53  can  be 
rewritten  as 

r'{t)  =  V2Paoco8(2Trfot  -b  ^0)  +  y/^aiRc{oc)cos{0i  —  6o)co8{2Trfot  -b  0o) 

—V^aiRc{a)8in{0i  —  9o)8in{2Tfot  +  6o)  (54) 

Applying  vector  analysis,  as  illustrated  in  Figure  17,  Equation  54  can  be  written  as  [6] 

/(<)  =  C{a)co8(2irfot  +  ^0  +  <i>e)  (55) 


where  C{a)  is  the  magnitude  of  the  combined  signals  and  4>e  is  the  phase  deviation  due  to 
the  reflected  signal  emd  is  given  by 


<t>e 


tan  ^ 
tan~^ 


V^CilRcioc)8in{9i  —  ^0) 
V^ciQ  +  V^o-iRc{oi)co8{9i  —  9o) 

—Rc{a)8in{9o  —  ^i)  1 

oo/ai  -b  Rcia)co8{9o  -  0i)J 


(56) 


where  0o  ~  ^1  =  2TfoC(Tc.  Prom  Equation  55,  the  phase  of  the  signal  from  the  PLL,  ^3, 
can  be  expressed  as  ^3  =  ^0  +  ^e- 

Equation  56  provides  insight  into  the  effects  of  the  the  code  phase  and  amplitude 
of  the  reflected  signal.  Figure  18  plots  (j>e  versus  a  for  oq/oi  =  2  euid  foTc  =  10.  The 
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Figure  17  Vector  representation  of  r'{t).  Note:  ip  =  9i  —  00- 

phase  error  envelope  shares  properties  with  the  triangular  code  auto  correlation  function, 
Rc,  of  Equation  8.  Also,  the  carrier  phase  of  the  PLL  generated  signed  is  in-phase  with  the 
direct-path  carrier  (i.e.,  (pe  =  0)  whenever  a  is  at  delays  that  corresponds  to  2'KfcaTc  =  Trn 
where  n  is  an  integer.  This  corresponds  with  the  reflected  signal  being  in-phase  or  180 
degrees  out  of  phase  with  the  direct-path  signal. 

3.3  MCTL  Operation 

Mixing  the  received  signal,  rmp{t),  with  the  signal  from  the  PLL,  2cos{2irfot  -I-  03), 
produces 

r(<)  =  V^Paoc(t  —  ro)cos(0o  —  03) V^2Paoc(t  —  ro)cos(47r/ot -I- 00  +  03) 

-|-V^aic(t  -  To  -  aTc)cos{9i  -  6z)  -h  V^Paic(t  —  To  —  aTc)cos{4irfot  -I-  0i  -I-  63) 
+2n{t)cos{2wfot  -h  03)  (57) 

which  enters  the  MCEU  and  the  MCTL.  The  ‘early-late’  gates  of  the  MCTL  operate  similar 
to  the  ‘early’  and  ‘late’  branches  of  a  coherent  DLL,  in  which  the  received  signal  is  mixed 
and  filtered  to  baseband.  eMRDLL  is  unique  in  that  it  implements  a  floating  ‘early-late’ 
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Figure  18  Local  carrier  phase  error,  ^e,  for  oo/oi  =  2  and  foTc  =  10. 

gate  code  tra.cking  loop  for  the  reflected  signal  in  addition  to  the  ‘early-late’  gate  tradsing 
loop  for  the  direct  path  signal.  The  floating  ‘early-late’  code  tracking  loop  uses  estimates 
of  the  multipath  parameters  provided  by  the  MCEU  to  synchronize  the  correlators  of  the 
‘early’  and  ‘late’  branches.  Unlike  the  traditional  DLL,  the  tracking  loops  of  the  MCTL 
are  able  to  synchronize  to  the  reflected  signal  as  well  as  the  direct  path  signal  to  reduce 
code  phase  tracking  error. 


‘early’-‘late’  gate 

Figure  19  ‘Early’  ‘late’  gate  delay  lock  loop  block  diagram. 
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Since  the  signal,  r{t),  is  low  pass  filtered  in  the  MCTL  as  well  as  the  MCEU,  it  can 
be  simplified  as 


r(t)  =  xoc(t  —  To)  +  a;ic(t  —  To  —  aTc)  +  n'{t) 


(58) 


where 


ajo  =  V2Paocos{0Q  —  6z) 

XI  =  ‘^/^a\cos{9i  —  9z)  (59) 

and  n'{t)  is  lowpass,  zero-mean,  AWGN  given  as 

n'ify  =  y/2[ni(t)cos(9z)  -l-nQ(t)sm(03)].  (60) 


Upon  entering  into  the  ‘early’  and  ‘late’  branches  of  the  tracking  loop,  the  signal,  r{t),  is 
mixed  with  the  locally  generated  spreading  code  defined  as 

CAit-ro-t3kTc)  =  cit-fo-^kTc-^)-c{t-fo-0kTc  +  ^)  (61) 

where  MRDLL  has  traditionally  implemented  a  one  chip  spacing,  A  =  1.0.  After  mixing, 
the  signals  in  each  branch  are  lowpass  filtered. 


3.3.1  Direct  Path  Channel.  For  high  processing  gains,  such  as  GPS,  and  small 
loop  bandwidths  {Bl  <10  Hz),  the  filtered  signals  of  the  ‘early’  and  ‘late’  branch  of  the 
direct  path  channel,  for  whidi  /3fcTc  =  0,  can  be  expressed  as  [9,11] 


yd,e{t)  = 


a:oc(t  -  ro)c(t  -  fo  -t-  yTc)  +  ®ic(t  -To-  aTc)c{t  -To  +  yTc) 


+  n'{t)c{t  —  To  +  -tTc)] 

^  JLPF 

y)  [(«  +  «  +  y) 


•—  XqRq 


n'i 


+  \  n  {t)c{t  -  To  +  jTc) 


LPF 


(62) 
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and 


yd,i{t)  = 


Xoc(t  —  T<,)c(t  -  To  —  ^Tc)  +  Xic{t  —  To  —  aTc)c{t  -To  —  -^Tc) 


+n\t)c{t  -To-  —Tc) 


=  XqRc 


iLPF 

+  [n'(t)c(t-fo-yrc)l 
L  ^  J  LPF 


\ 


respectively.  Differencing  the  ‘early’  and  ‘late’  branches,  yd,i  —  yd,e^  produces 

yd{t)  =  ®o{iJc  [(«-!) To  -Rc  (^  +  1)^'=]} 

+®1  •|iic  +  o:  —  Tc  —  Rc  +  a  +  Tcj  | 
+  [n'(t)cAit-To)]j^pp 


(63) 


(64) 


After  the  ‘early-late’  differencing  operation,  the  signal  enters  the  gain-phase  correlators 
upon  which  the  signal  is  mixed  with  the  MCEU  signal  parameter  estimate  of  ®o;  the 
output  of  this  direct-path  arm  can  be  expressed  eis 

Zd{ti  S)  =  ®o®o  l-Rc  -  -Rc  Tcj  I 

-I-S0®!  I-Rc  +  a  -  Y^  Tcj  -  Rc  -b  a  +  y^  Tc  | 

-|-®0  [n'(t)cA(t  -  To)]  j^pp .  (65) 
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3.3.2  Multipath  Channel.  Similarly,  the  mixed  signals  for  the  ‘early’  and  ‘late’ 
branch  of  the  reflected  path  channel  can  be  expressed  as  follows 

yT,e{t)  =  ®oc(i  -  To)c(i  -  To  - /JfeTc  +  yTc) 

A 

+  Xic{t  -  To  -  aTc)c{t  -  To  -  /dfeTc  +  —Tc) 

+  n'{t)c{t  —  To  —  l3kTc  +  -rlc)] 

^  JLPF 

=  xoRo  +  + 

+  n'{t)c{t  -  To  -  PkTc  +  4rc)]  (66) 

^  J  LPF 


A 

yrj{t)  =  Xoc{t  -  To)c{t  -  To  -  PkTc  -  —Tc) 

+xic{t  -  To  -  aTc)c{t  -To-  PkTc  -  —Tc) 

+  n'(t)c(t  -  To  -  pkTc  -  ^Tc)] 

^  JLPF 

=  XqRc  Oi 

+  n'(t)c(t  -  To  -  ^kTc  -  “Tc)  I 

^  J  LPF 


respectively.  Differencing  the  ‘early’  and  ‘late’  branches,  yr^i  —yr,ei  produces 
2/7.  =  ajo  ^-Rc  0£  —  O'  +  ^cj  ^ 

+»i{iee[(«-f)r,]-it.[(i  +  |)r,]} 

+  [n'(t)cA(t  -  To  -  /3kTc)]j^pp .  (68) 

After  the  ‘early-late’  differencing  operation,  the  signal  enters  the  gain-phase  correlators 
upon  which  the  signal  is  mixed  with  the  MCEU  signal  parameter  estimate  ®i.  The  output 


38 


of  the  multipath  arm  can  be  expressed  as 


-iie  («  +  f)r.]| 

+®l  [n'{t)cAit  -To-  /3kTc)]  (69) 


-Ro  ^6-a  +  l^re  I 


S.3.3  Discriminator  Output.  Under  the  assumption  that  the  MCEU  provides 
perfect  estimates  of  sq  and  »i  to  the  gain-phase  correlators,  the  respective  outputs  of  the 
direct-path  and  multipath  branches  are  [6] 

Zd{t,  6)  =  XqD{8)  +  xoxiD{S  -I-  a)  -I-  »o  [n'(t)cA(t  -  To)] 

Zr{t,  6)  =  xoxiD{6  -  a)  +  XiD{S)  +  xi  [n'(t)cA(t  -  h  -  ^kTc)]ipp  (70) 

where  the  MCTL  D  —  curve  is  defined  as 

i>(5)^iZc[(5-|)ro  (71) 

Summing  the  outputs  of  the  direct-path  channel  and  the  multipath  channel  produces  the 
code  tracking  discriminator/error  signal 

e(t)  =  (®o  +  ®i)7?(0  +  ®0®l[7?(^  +  a) -H  iP(5  —  a)] 

-|-»o  [n'(t)cA(t  -  To)]  ^pp  +  xi  [n'(t)cA(t  -  To  -  ^kTc)] p  (72) 

which  can  be  expressed  as  [6] 

e(t,5)  =  5(5)-|-ne(t)  (73) 

where  S{6)  represents  the  MRDLL  S-curve  given  as 

S{6)  =  (®o  +  Xi)D{S)  -h  X(iXi[D{S  -h  a)  +  D{S  -  a)]  (74) 
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and  the  discriminator  output  noise,  ne{t),  is  given  as 


ne{t)  =  XQ  [n'(t)cA(t  -  To)]j^pp  +  XI  [n'(t)cA(t  -  To  -  PkTc)]j^pp  (75) 

Figure  20  provides  a  plot  of  the  discriminator  components  of  the  direct-path  channel 
and  the  multipath  channel.  The  figure  also  shows  that  the  sum  of  the  two  channels  produces 
the  discriminator  output  which  has  an  error  signal  equal  to  zero,  under  the  condition  of 
perfect  estimation. 


Figure  20  Typical  MRDLL  S-curve  discriminator, 

3.3.4  Determination  of  Noise,  ne(t).  The  noise  components  of  ne{t)  were  ob¬ 
tained  under  the  assumption  that  the  LPFs  in  each  of  the  ‘early-late’  gates  were  ideal. 
Alternatively,  the  noise  component,  ne(t),  can  be  represented  as 

ne(t)  =  xond{t)  +  xiUrit)  (76) 

where  nd{t)  and  nr{t)  are  described  by  [6] 

njb(t)  =  ^/2  |^n{(t)cos(^3)  +  (t)sm(d3)j 

”*(<)  =  ^kTc)*hipF{t) 

=  nQ{t)cAit  -  To  -  l3kTc)  *  hLPF{t)  (77) 
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where  k  =  d  oi  r,*  denotes  convolution,  and  hippit)  is  the  impulse  response  of  the  LPFs. 
For  large  code  period,  IV  1,  and  small  loop  bandwidth,  Bi  <  10  fl’z,  the  amplitude  of 
the  two-sided  PSD,  Sn'k^  of  n{(t)  and  n^{t)  is  approximately  Nq  in  units  W/Hz  [9].  As 
such,  the  output  of  the  noise  of  the  k*^  branch  has  a  two-sided  PSD 


Sn,U)  = 


{ 


21'ro  \f\<LPF 
0  elsewhere 


(78) 


Under  the  condition  that  the  MCEU  provides  perfect  estimates  to  the  gain-phase  correla¬ 
tors,  the  combined  PSD  noise  term  of  Zi(t)  and  Zr{t)  is 


SnSf) 


I 


2iVo(®o  +  ®i) 

0 


I/I  <  LPF 

elsewhere 


(79) 


3.3.5  eMRDLL  Linear  Model  and  Tracking  Performance.  The  operation  of  the 
VCC  within  the  MCTL  provides  the  same  performance  characteristics  as  described  in 
Chapter  2  for  the  NCDLL.  Therefore,  the  output  phase  of  the  VCC  can  be  expressed  as 


To 


Ko  f  e{a,6)*fi\)dX 
Jo 

Ko  f  f  e{a,6)f{a  —  X)d\da. 
Jo  J — oo 


(80) 


where  the  output  of  the  loop  filter  is  described  by  its  impulse  response,  /(t),  is  convolved 
with  the  input  signal,  e{t,6).  Substituting  Equation  73  into  Equation  80  provides  the 
equation  for  the  nonlinear  model 

tM  =  Ko  f*  r  {S{6)  +  ne(X))  fia  -  X)dXda  (81) 

-‘c  Jo  J—oo 


where  6  =  (ro(t)  —  fo{t))/Tc. 

For  small  tracking  errors,  the  linear  region  of  the  S-curve  can  be  expressed  as  [6] 

5(5)  =  A6  =  A  (82) 
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where  A  is  defined  as  the  slope  of  the  lineeur  operating  region  for  small  S.  Substituting 
Equation  82  into  the  nonlinear  equation,  Equation  81,  results  in 


^o(<) 

Tc 


=  AKc, 


r*  &  Y To(a)  -  -ro(a) 

i/o  </— (X  ^ 


) 


+ 


neja) 

A 


f{\-a)dadX 


(83) 


which,  in  the  Laplace  domain,  can  be  represented  as 

fojs)  _  Tojs)  -  fo(8) 

Tc  Ta 


AK, 


F{8) 


(84) 


where  F{s)  is  the  Laplace  transform  of  f(t).  Therefore,  the  linear  equivalent  circuit  model 
for  the  MCTL  can  be  represented  by  replacing  the  terms  and  n«(t)  from  the  model  in 
Figure  6  with  A  and  ne(t)  as  shown  in  Figure  21. 


Figure  21  eMRDLL  linear  equivalent  circuit. 


Solving  for  To{s)/to{s)  given  Equation  84  results  in  the  closed  loop  transfer  function 


H(s\  4 

^  ^  To(s)  s  +  AKoFis) 


(85) 


which  characterizes  the  eMRDLL  for  fixed  estimates  d,  »o>  aJid  *!•  Choosing  an  active  lead- 
lag  loop  filter,  the  transfer  function  of  Equation  85  can  be  transformed  into  the  classical 
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model  representation  given  as  [2, 9] 


H{8) 


2CuJnS  + 
s2  +  2C(t>nS  +  Ul 


(86) 


where  the  loop  natural  frequency,  a;„  in  radfsec,  and  the  (imitless)  damping  factor,  axe 
defined  as 


Wn  = 


(87) 


and 


C  =  yWn  (88) 

respectively.  As  with  the  NCDLL,  the  single-sided  noise  equivalent  bandwidth,  Bi  in  Hz^ 
for  a  second-order  loop  with  an  active  lead-lag  loop  filter  can  be  described  as  [2,8] 

poo 

Bl  =  J  \H{j2irf)fdf 

=  !^(c+i).  m 


The  power  spectrum  of  the  tracking  jitter  can  be  expressed  as  [2,9] 


Ssif)  =  \HU2Tf)fSne'if)  (90) 

where  H{s)  is  the  closed-loop  transfer  function  defined  in  Equation  86,  and  Sne'if)  = 
is  the  two-sided  power  spectrum  of  the  Gaussian  noise  process  at  the  input  to 
the  loop  model.  Ss{f)  is  approximately  fiat  over  the  loop  bandwidth,  Bl,  so  that  the 
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variance  of  6  can  be  expressed  as 


/oo 

S6if)df 

‘OO 

=  ^Swy)|ffO'2>r/)|V 
4(go  +  icj)^o„ 

”  a7 


(91) 


where  the  integral  on  the  right  side  of  the  equation  is  defined  as  the  two-sided  noise 
bandwidth  Bl  in  Hz. 

As  stated  in  Chapter  2,  the  natural  frequency,  Un,  and  the  damping  factor,  C,  de¬ 
termine  the  tracking  performance  of  the  code  tracking  loop  and  the  noise  loop  bandwidth, 
Bl,  determines  the  amount  of  noise  power  that  affects  the  tracking  loop.  As  shown  in 
Equations  87  and  88,  the  slope  of  the  lineeir  region  of  the  S-curve,  A,  which  varies  for 
different  values  of  amplitude  eind  delay  of  the  received  multipath  signal,  sdso  determines 
the  response  characteristics  of  the  loop  via  and  C-  As  such,  the  changing  value  of  A,  as 
tabulated  in  Table  2,  will  also  impact  tracking  performance. 


3.4  ALC  Operation 

Because  the  transient  response  of  the  loop  is  dependent  upon  the  parameters  of  the 
reflected  signal,  a  variable  gain  is  introduced  in  an  attempt  to  eliminate  the  dependency. 
This  is  accomplished  by  preceding  the  loop  filter  with  a  variable  gain,  as  done  in  the  work 
of  Laxton  and  DeVilbiss  [6,7]. 

As  shown  in  Figure  22,  the  ALC  is  a  second  order  loop  controller  preceded  by  a 
variable  gain.  The  variable  gain  allows  the  designer  to  fix  the  dynamic  response  of  the 
MCTL  within  a  linear  operating  region.  The  ALC  is  a  unique  element  not  present  in  the 
RDLL  design  of  [10].  The  ALC  uses  the  parameters  provided  by  the  MCEU  to  control 
tracking  performance  relative  to  the  discriminator  output  [6,7]. 
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a  xi  ajQ 


■►(to  VCC) 


I _ j 

ALC 

Figure  22  Adaptive  loop  controller  (ALC)  block  diagram. 

Assuming  a  2"**  order  active  lead-lag  loop  filter,  given  as 

Fa{s)  =  (92) 

8 

the  filter  shall  be  preceded  by  a  variable  gain,  Ka  —  'i-fri.  This  provides  a  composite  loop 
transfer  function  of 


1 

1 

1 

1 

II 

.31- 

Fa  =  ^ 

1  ^ 

I 

1 

1 

Adaptive  Gain 

Loop  Filter 

FW  =  jr^2£±i 

8 

_  r2S  +  l 
ri8 


(93) 


As  shown  in  Table  2  and  illustrated  in  Figures  23  and  24,  the  slope  of  the  linear  operating 
region.  A,  varies  with  the  signal  parameters  of  rjnp{t)-  The  variable  gain  of  the  ALC  allows 
adaptation  to  changing  signal  parameters  so  as  to  maintain  a  fixed  natural  firequency, 
and,  therefore,  a  fixed  damping  ratio  and  noise  equivalent  loop  bandwidth,  Bi. 


The  ALC  uses  received  signal  estimates  d,  xq,  and  ®i  firom  the  MCEU  to  determine 
an  estimate  of  the  slope  of  the  linear  tracking  region.  A,  based  on  the  equations  provided 
in  Table  2  [6,7].  The  estimate.  A,  is  then  used  to  adjust  the  gain 


Ka  = 


1 

n  ~  AKo 


(94) 


where  a;„  and  Ko  are  fixed  parameters.  The  coefficient,  T2,  in  the  loop  filter  remains  fixed 
based  on  Equation  88  which  is  specified  by  the  designer  based  on  desired  performance 
characteristics. 
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Table  2  MCTL  linear  operating  region. 


Gain  =  (A/2) 

Multipath  Delay  Range 

Tracking  Error  Range 

»o  +  ®i  +  2»o®i 

0<0.5 

-0.5  +  o  <  5  <  0.5  -  a 

®o  +  ®1  +  5®0®1 

a  =  0.5 

-0.5  <  5  <  0.5 

®0  +  ®l  -  ®0®1 

0.5  <  1.0 

0.5  -  a  <  5  <  -0.5  +  o 

XQ  +  xf-  ®0®1 

1.0  <  1.5 

-1.5  +  0  <5<  1.5-0 

®0  +  ®l  “  5®0®1 

0  =  1.5 

-0.5  <  5  <  0.5 

Figure  23  eMRDLL  S-curve  with  xq  =  1,  ®i  =  0.5,  and  a  =  0.5. 

5.5  MCEU  Operation 

The  MCEU  consists  of  a  bank  of  M  correlators  followed  by  a  ML  estimation  block 
as  shown  in  Figure  25.  Each  correlator  consists  of  a  mixer  followed  by  a  low  pass  filter, 
which  acts  as  an  integrator  when  high  processing  gains  are  used.  The  estimation  block 
generates  estimates  of  the  reflected  signal  delay  relative  to  the  direct  path  signal,  as  well  as 
coefficients  of  the  LOS  and  reflected  baseband  signals;  these  estimates  Eire  fed  to  the  MCTL 
and  the  ALC.  The  MCEU  is  capable  of  estimating  the  parameters  of  multiple  reflected 
signals.  As  such,  the  received  signal,  r(t),  entering  the  MCEU  is  represented  as 

N 

r(t)  =  xoc{t  —  To)  +  ^  ®jc(t  -To  -  ocxTc)  +  n(t)  (95) 

t=i 
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Figure  24  MRDLL  S-curve  with  xq  =  1,  xi  =  0.5,  and  a  =  0.7. 

where  ®i  =  \/2Po<cos(^<  —  Op),  emd  N  represents  the  number  of  reflected  signeds.^ 

The  received  signal  enters  the  bank  of  correlators,  where  each  arm  is  uniformly 
spaced,  at  through  the  range  of  [0,1.5].  Under  perfect  tracking  conditions,  where 
fo  =  To,  the  sampled  output  of  the  correlator  can  be  expressed  as 

N 

Rk{nTs)  =  xoRc{/3k)  +  ^  XiRciai  -  f3k)  (96) 

t=i 

where  j5k  is  the  delay  of  the  respective  correlator  arm,  relative  to  f©,  and  represents  the 
sampling  period  of  the  estimator. 

Let  a  be  defined  as 


a  =  [ai  0(2  ’  • 


(97) 


The  sampled  outputs  from  the  bank  of  M  correlators  entering  the  estimator  block  are 
collectively 


R  =  [iZo  i2i  •  •  •  Rm-\Y 


(98) 


^dp  represents  the  phase  of  the  signal  generated  from  the  PLL,  which  was  previously  represented  as  ^3. 
It  has  been  changed  for  this  section  to  remove  ambiguity  with  the  number  of  multipath  signals. 
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MCEU 


Figure  25  Multiple-correlator  estimation  unit  (MCEU)  block  diagram. 

which  can  be  expressed  as  the  hnear  statistical  data  model,  R  =  H(a)x  -1-  v.  H(a)  is  an 
M  y.  N  +  1  observation/regressor  matrix  where  M  is  the  number  of  correlator  arms  and 
iV  -h  1  is  the  total  number  of  signal  paths.  H(a)  can  be  represented  as 

RcWo)  Rcioii  —  Pq)  •••  Rci^N  -  M 
H(a)=  :  i  i  (99) 

_  Rc{Pm-i)  Rc(«1  -  Pm-i)  •  •  •  Rc(oin  —  Pm-i)  _ 
and  the  composite  signals  of  the  correlator  branches  are  collectively 

X  =  [»o  •  •  •  ®  jv]^*  (100) 

The  observation  noise  vector,  v,  is  modeled  [5]  as  zero  mean  Gaussian  with  covariance 
Cv.  To  determine  an  expression  for  Cy,  the  spreading  code  is  considered  a  random  binary 
process^,  independent  of  the  observation  noise.  The  continuous-time  cross  correlation 
between  the  and  correlator  noise  outputs  is 
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which,  for  nctrrowband  low  pass  filtered  noise,  can  be  approximated  as  [6] 

Rij  (r)  «  i2c(T  +  A/3ij)<Tl  (102) 

where  A/3jj  =  —  Pj  and  is  the  variance  of  each  of  the  correlator  noise  outputs. 

As  such,  the  instantaneous  observation  noise  covariance  matrix,  C-„  =  E  [v(t)v'^(t)]  — 
E[vit)]E[y'^it)]  ,  can  be  expressed  as 


Cv  =  cr?C 


(103) 


where  C  is  the  M  x  M  symmetric,  Toeplitz  correlation  matrix  [6] 


C  = 


1  i?c(/3o-/5i) 

Rc{l3i-/3o)  1 

Rc(J^m-i  —  M 


•  •  •  Rc{/3o  -  ^M-l) 

RcifiM-2-  Pm-i) 
RcifiM-1  -  Pm-2)  1 


•  (104) 


For  a  general  linear  data  model,  R  =  H(a)x  +  v,  where  H(a)  is  an  AT  x  p  matrix  of 
known  structure  and  v  ~  Ar(0,  Cv)  is  a  noise  vector  of  dimension  AT  x  1,  the  mccisurement 
probability  distribution  function  (PDF)  can  be  expressed  as  [5] 


p(R)  = 


(27r)^/2|Cv|l/2 


exp' 


(-l/2(R-H(a)x)'rCv-‘(H'-H(a)x)] 


(105) 


The  maximum  likelihood  estimate  of  x,  xml>  is  found  by  minimizing  the  quadratic  func¬ 
tion  [5] 


F(®)  =  (R  -  H(a)x)^Cv"HR  -  H(a)x). 


(106) 


The  first  order  necessary  condition  to  optimize  F(a)  is  that  dV{x)/dx  =  0.  For 
symmetric 


gF(a) 

dx 


=  2H(a)^Cv"^(R  -  H(a)x) 


(107) 
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Setting  dV{x)ldx  equal  to  zero  yields  (independent  of 

H(af  C-^il  =  H(af  C-lH(a)®ML  (108) 

If  H(a)  is  of  full  rank  and  N  <  M  —  1,  then 

XML  =  (H(a)C-^H(a))-^H(a)^C-iR  (109) 

Because  C“^  is  positive  definite,  the  existence  of  (H(a)C“'H(a))“^  is  sufficient  to  guar¬ 
antee  that  Xml  is  a  global  (unconstrained)  minimum.  The  MLE,  xml>  has  been  shown 
to  be  the  minimum  variance  unbiased  (MVU)  estimator  as  well  as  an  efficient  estimator 
by  [5].  The  PDF  of  xml  is  given  as  iV(x,  (H(a)^Cv“^H(a))“^). 

Xml,  via  H(a),  is  actually  a  function  of  a,-  G  (0, 1.5],  which  does  introduce  constraints 
to  the  problem.  If  the  primary  path  signal  is  perfectly  synchronized  with  the  receiver 
generated  replica,  then  a,-  <  0  is  not  physically  possible  and  a*  >  1.5  corresponds  with  a 
delayed  signal  which  heis  no  impact  upon  the  measurements,  R  [19]. 

An  optimization  to  find  the  minimizing  set  of  Oi  €  (0,1.5]  can  be  obtained  by  sub¬ 
stituting  Xml  into  Equation  106  as  follows  (Note:  for  >  0,  minimizing  F(q:)  s  (7^V{a) 
is  equivalent  to  minimizing  V{a).) 

Via)  =  (R-Hx^)C-i(R-Hx)|,=i^,(„) 

=  (R^  -  R^C-1H(H^C-1H)-1H^)C-HR  -  H(H^C-iH)-iH^C-iR) 

=  R^C-^R  -  R^C-^H(H^C-1H)-^H^C-^R 

=  (C-^R)^[C-H(H^C-^H)-^H^](C-iR).  (110) 

C  and  R  are  known  and  measured,  respectively,  and  the  vector  of  unknowns  is 
a  =  [ai  aa^’Q^iv]^  of  Equation  97,  which  is  needed  to  specify  H(a)  of  Equation  113. 
A  minimization  of  V’(a)  yields  ocml  which,  in  turn,  determines  xml  via  Equation  109. 
Common  numerical  optimization  techniques  for  minimizing  Equation  110  are  quasi-Newton 
methods  such  as  sequential  quadratic  programming  (SQP).  There  are  risks  associated  with 
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these  iterative  techniques:  the  iteration  may  not  converge,  the  iteration  may  converge  to 
a  local  minimum  as  opposed  to  the  global  minimum,  or  the  signal-to-noise  ratio  may  be 
insufficient  for  meaningful  results. 

3.5.1  Analysis  of  MCEU  for  a  Single  Reflection.  Returning  to  the  case  of  one 
reflected  signal,  the  signal  entering  the  MCEU  can  be  represented  eis 

r{t)  =  xoc{t  —  To)  +  »ic(<  —  To  —  aTc)  +  n'{t)  (111) 

where  xq  =  y/^aocos(6o  —  63)  and  xi  =  \/^oicos(^i  —  63).  Let  Vk(t)  represent  the  noise 
of  the  MCEU  correlator.  Assuming  perfect  code  tracking,  the  correlator  outputs  are 

Rk{t)  =  xoRciflk)  +  »li2c(ai  -  flk)  +  Vk(t)  (112) 


for  k  =  0, 1, . . .  ,  M  —  1  which  can  be  represented  collectively  by  the  lineeu:  data  model 
R  =  H(a)x  +  V,  where  x  =  [»o  ®l]^  (Equation  100  for  N=l), 


H(a)  = 


RcWo) 

i2c(Q!l  -  /?o) 

RcWi) 

Rcioci  -  fll) 

Rc{fiM-\) 

Rcioci-  Am-i) 

(113) 


R  is  given  in  Equation  98,  and  a  is  given  in  Equation  97  for  JV  =  1.  When  fo  =  To,  the 
flrst  column  of  H(a)  represents  the  cross  correlation  with  the  direct  signal  and  the  second 
colunm  represents  the  cross  correlation  with  reflected  signal. 


3.5. 1.1  Numerical  Optimization.  Given  that  the  samples,  R,  obtained  at 
the  output  of  the  bank  of  correlators  conform  to  the  statistics  of  Equation  105,  the  cost 
function  V  (a)  of  Equation  110  can  be  minimized  numerically.  Quasi-Newton  methods  build 
up  curvature  information  of  V  (a)  at  each  iteration  to  formulate  a  quadratic  model  problem 
(provided  in  Matlab  Optimization  Toolbox)  from  which  the  optimal  solution  is  determined 
when  the  partial  derivatives  of  the  quadratic  function  go  to  zero.  By  using  the  observed 
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behavior  of  the  quadratic  function  to  build  up  curvature  information,  an  approximation 
of  the  minimum  can  be  determined  with  much  fewer  operations  than  traditional  Newton- 
type  methods  which  proceed  in  the  direction  of  descent  of  the  quaxiratic  function  using  a 
line  search  method  to  locate  the  minimum  after  several  iterations.  Sequential  Quadratic 
Programming  (SQP)  minimizes  Newton’s  method  for  constrained  optimization  using  a 
quasi-Newton  updating  method  to  approximate  the  Lagrangian  function  of  V’(a). 

3. 5. 1.2  Theoretical  Analysis  ofV(a)  for  a  Single  Reflection.  The  success 
of  any  SQP  optimization  routine  depends  on  the  algorithm’s  ability  to  determine  the 
curvature  of  the  given  quadratic  function  from  which  it  can  converge  to  the  minimum. 
Convergence  to  global  minimum  of  a  quadratic  function  can  become  very  difl&cult  given  a 
series  of  local  minimums.  Figure  26  provides  1^(0:)  versus  all  possible  values  of  a  €  [0,1.5] 
for  a  given  multipath  signal  with  delay  a.  As  shown  in  the  figure,  the  absolute  or  global 
minimum  is  indeed  at  the  location  of  the  multipath  delay.  However,  the  curvature  of  V (a) 
varies  through  the  range  of  as.  The  changes  of  slope  (i.e.,  the  ridges)  of  V(a)  occur  at 
the  location  of  the  arms  of  the  correlators,  flkTc.  In  Figure  26  a  bank  of  M=ll  correlators 
spaced  evenly  at  increments  of  0.15  through  out  the  range  of  a.  In  addition  the  change 
of  slope  that  occurs  at  0.05,  it  is  seen  that  there  are  M  —  1  curvature  changes  that  occur 
at  the  delays  corresponding  to  correlator  locations  interior  to  (0,1.5).  Figure  27  provides 
V(a)  versus  a  €  (0,1.5)  for  M=16.  Once  again,  there  M  —  1  changes  of  curvature,  at  the 
location  of  the  interior  correlator  arms. 

Because  the  slope  of  function  is  not  defined  at  the  location  of  the  ridges,  convergence 
problems  can  occur,  as  previously  described.  An  SQP  algorithm  could  potentially  estimate 
the  curvature  of  the  slope  of  V  (a)  incorrectly  which  would  corrupt  the  estimate  of  the  delay, 
d  of  the  reflected  signal.  As  such,  the  correlator  spacing,  flkTc,  and  the  search  algorithm 
employed  should  be  considered  carefully. 

A  previous  design  incorporated  correlators  spaced  apart  by  O.lTc  through  the  range 
of  possible  a  G  [0, 1.5],  which  resulted  in  a  bank  of  16  correlators  [6].  We  have  empirically 
determined  that  implementing  a  bank  of  11  correlators,  evenly  spaced  apart  by  0.15rc, 
improves  the  performance  of  Matlab’s  constr.m  SQP  algorithm.  By  reducing  the  number 
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Plot  of  V(alpha)  values  (M-11) 


Figure  26  V’(a)  given  a  =  0.7  where  M=ll  correlators  Eire  spaced  at  increments  of  0.15 

through  the  range  [0,1-5]. 

of  correlators,  fewer  changes  of  slope  of  the  function,  V^(q;),  occur  which  improves  the 
SQP  algorithm’s  ability  to  converge;  this  simultaneously  reduces  the  amoimt  of  hardware 
required.  It  was  found  that  reducing  the  number  of  correlators  further,  led  to  degraded 
SQP  algorithm  performance. 

Convergence  problems  have  also  been  determined  to  occur  for  reflected  signals  with 
delays  very  close  to  the  direct  path  signal.  As  the  reflected  signal  becomes  very  close  to 
the  direct  path  signal,  the  columns  of  H(a)  become  increasingly  linearly  dependent  and 
H(a)  becomes  rank  deflcient.  As  such,  it  is  more  difiicult  to  numerically  estimate  the 
delays  of  reflected  signals  with  delays  very  close  to  the  direct  path  signal.  Therefore,  the 
SQP  constraint  function  implemented  was  limited  to  estimating  reflected  signal  delays  of 
a  6  [0.05,1.5].  Figure  28  provides  V’(a)  versus  a  where  the  multipath  delay  is  0.05.  The 
global  minimum  of  the  function  is  located  at  0.05  which  constitutes  the  edge  of  the  search 
region.  However,  if  allowed  to  seeirch  the  region  [0,0.05]  the  global  minimum  may  not  be 
located  because  the  slope  through  this  region  is  essentially  zero.  Figure  29  provides  V (a) 
versus  a  where  the  multipath  delay  is  0.1.  Once  again,  the  global  minimum  of  the  function 
is  located  at  the  delay  of  the  multipath  signal. 
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Plot  of  V{elptia)  values  (M-1 6) 


Figure  27  V^(a)  given  a  =  0.7  where  M=16  correlators  eure  spaced  at  increments  of  0.1 

through  the  range  [0,1-5]. 

Figures  30  and  31  provide  V(a)  versus  a  for  the  region  where  the  delays  of  the 
multipath  signal  are  1.5  and  1.4,  respectively.  As  shown  in  the  figures,  convergence  to  the 
global  minimum  does  not  constitute  a  problem  for  these  delays. 

3. 5. 1.3  Search  method  for  global  minimum  ofV{a)  for  single  reflection.  In 
order  to  reduce  convergence  problems,  a  partitioned  (ID)  search  technique  can  be  used 
to  find  the  global  minimum  of  F(a).  Since  the  undefined  points/ridges  are  known  to 
be  located  at  the  placement  of  the  correlator  arms,  a  search  for  the  local  minimum  can 
be  performed  for  the  intervals  between  each  of  the  arms.  The  local  minimum  for  each 
interval  can  be  determined  by  some  numerical  optimization  routine  and  then  compared 
to  determine  the  global  minimum.  This  method  requires  that  the  optimization  routine  be 
performed  for  each  interval;  however,  convergence  within  each  interval  will  require  fewer 
iterations  than  a  search  across  the  entire  remge  of  a,  since  the  quadratic  function  through 
each  interval  is  smooth  (with  no  undefined  points/ridges).  The  global  minimum  could 
possibly  be  converged  on  more  quickly  by  searching  each  interved  using  parallel  processing 
such  that  the  local  minimums  through  the  peurtitions  of  V  (a)  are  determined  concurrently. 
This  approach  can  be  generalized  for  more  than  a  single  reflection. 

Table  3  provides  the  results  of  the  search  algorithm  employed  using  the  constr.m 
fimction  in  Matlab  for  a  single  multipath  delay  and  M=ll  correlators.  The  results  provided 
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PkitotV(alpha)valuM(M-11) 


Figure  28  V (a)  given  a  =  0.05  where  M=ll  correlators  are  spaced  at  increments  of  0.15 

through  the  range  [0,1'5].. 

in  the  table,  for  a  in  uniformly  spaced  increments  of  0.05,  indicate  that  the  algorithm  was 
able  to  estimate  each  of  the  signal  parameters  with  an  accuracy  >  99.95%. 

3.5.2  Analysis  of  MCE  U  for  Multiple  Reflections.  For  the  case  of  two  reflected 
signals,  x  =  [so  si  *2]^  and  H(a)  is  an  M  x  3  matrix  with  the  first  column  representing 
the  cross  correlation  with  the  direct  signal  and  the  second  emd  third  columns  representing 
the  cross  correlation  with  each  of  the  reflected  signals.  Figure  32  provides  V’(a)  versus 
all  possible  values  of  ai  and  0:2  through  the  range  ai,Q!2  €  [0, 1.5].  The  ridges  or  changes 
of  curvature  are  apparent  along  the  direction  of  ai  and  02.  The  figure,  as  expected,  has 
symmetry  about  the  diagoncil  of  ai  =  02]  the  elements  of  V’(a)  are  undefined  at  these 
locations.®  The  physical  interpretation  of  ai  =  0:2  is  that  only  a  single  reflection  exists. 

Figure  33  represents  a  slice  taken  at  the  known  location,  a  =  0.2,  of  one  of  the 
given  multipath  delays.  This  figure  exhibits  the  same  characteristics  as  the  case  for  one 
reflection.^  Ridges  occur  at  the  locations  of  the  arms,  pkTa  as  previously  described. 
In  addition,  the  absolute  or  global  minimum  is  located  at  the  point  of  the  other  given 
multipath  delay. 

^The  diagonal  elements  of  1^(0:)  have  been  set  to  some  arbitrary  value  in  the  2D  figures  presented 
although  they  are  undefined  because  H(a)  is  rank  deficient  when  ai  =02. 

'*The  point  at  0.2  should  be  ignored  since  it  as  the  location  of  the  diagonal  element. 
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0.036] 


Plot  Of  V(alpha)  values  (M-11) 


Figure  29  V (a)  given  a  =  0.1  where  M=ll  correlators  are  spaced  at  increments  of  0.15 

through  the  range  [0,1.5].. 

Figure  34  provides  V(a)  versus  Q'i,a2  G  (0,1.5]  given  two  multipath  delays  which 
differ  from  those  used  in  Figure  32.  Characteristics  similar  to  those  of  Figure  32,  such  as 
the  ridges,  can  be  recognized.  Figure  35  provides  a  close-up  of  the  location  of  the  given 
multipath  delays  of  Figure  34.  As  shown  in  the  figure,  the  global  minimum  is  indeed 
located  at  ai  =  0.5  and  aa  =  0.9.  Due  to  symmetry  of  F(q:)  about  the  diagonal,  a  like 
minimum  exists  at  ai  =  0.9  and  a2  =  0.5. 

Figure  36  provides  V(a)  versus  ai,Q'2  €  (0,1.5]  given  two  multipath  signals  with 
closely  spaced  delays,  ai  =  0.9  and  a2  =  1.0.  As  shown  in  the  figure,  V(a)  becomes 
very  fiat  neeur  the  location  of  the  multipath  delays,  which  may  cause  difficulties  for  an 
SQP  algorithm  to  converge  to  the  global  minimum.  In  fact,  without  a  search  routine,  the 
SQP  algorithm  implemented  in  Matlab  performed  poorly  with  multipath  delays  spaced 
relatively  close  together.  As  such,  a  search  technique  may  provide  increased  performance. 


3.5.2. 1  Search  method  for  global  minimum  ofV{a)  for  multiple  reflections. 

The  search  technique  applied  for  a  single  reflection  may  be  implemented  for  multiple 
reflections.  Implementing  a  search  technique  for  two  reflected  signals  requires  a  grid  (2D) 
search  technique  to  find  the  global  minimum.  Since  V (or)  is  symmetric  about  the  diagonal, 
only  grid  regions  on  one  side  of  the  diagonal  need  to  be  searched. 
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Plot  of  V(alpha)  values  (M-11) 


Figure  30  V^(a)  given  a  =  1.5  where  M=ll  correlators  are  spaced  at  increments  of  0.15 

through  the  range  [0,1.5]. 

This  symmetry  holds  true  for  more  than  two  reflected  signals  ets  well,  but  requires 
signiflcantly  increased  computational  burden  as  the  number  of  multipath  signals  increeises. 
The  number  of  search  regions  required  for  a  given  set  of  multipath  signals  can  be  expressed 
as  a  combination 


+  A  “  1  \ 

J 


where  k  equals  the  number  of  multipath  dimensions  (i.e.,  k-1  reflected  signals),  n  equals 
the  number  of  search  intervals  (in  one  dimension)  between  correlator  arms,  and  n>k. 

For  two  reflected  signals  with  M=ll  correlators,  the  number  of  search  regions  is 


(115) 


These  regions  were  searched  to  locate  the  global  minimum  for  the  given  set  of  multi- 
path  delays,  ai  and  0!2.  In  order  to  test  the  performance  of  the  SQP  search  algorithm  the 
true  ai  EUid  0:2  were  varied  by  increments  of  0.1  through  the  range  [.1,1.5]  which  provided 
105  test  cases.  Of  the  105  cases,  88  or  83.8%  of  the  cases  resulted  in  estimation  accuracy 
of  >  99.9%  for  the  psurameters  di,  0:2,  ®o>  ®l>  and  ®2*  The  results  of  the  remaining  17 
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Plot  of  V(alpha)  values  (M-11) 


Figure  31  V’(a)  given  a  =  1.4  where  M=ll  correlators  are  spaced  at  increments  of  0.15 

through  the  range  [0,1-5]. 

test  cases  are  provided  in  Table  4.  All  of  the  multipath  delays  presented  in  Table  4  axe 
spaced  <  0.2  apart  with  the  worst  results  occurring  for  delays  spaced  =  0.1.  Even  so,  the 
estimates  of  the  delays  were  within  0.1  chips  of  the  true  delay.  It  should  also  be  noted  that 
for  six  of  the  cases,  estimates  of  xi  and  ®2  were  very  high,  >  10.  This  information  could 
be  used  in  an  algorithm  to  indicate  very  poor  estimates. 
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Table  3  Table  of  trae  values  of  a,  xq,  and  ®i  with  corresponding  estimates 


a 

XQ 

XI 

a 

Xq 

XI 

0.0500 

1.0000 

-0.5000 

0.0500 

1.0000 

-0.5000 

0.1000 

1.0000 

0.5000 

0.1000 

1.0000 

0.5000 

0.1500 

1.0000 

-0.5000 

0.1500 

1.0000 

-0.5000 

0.2000 

1.0000 

0.5000 

0.2000 

1.0000 

0.5000 

0.2500 

1.0000 

-0.5000 

0.2500 

1.0000 

-0.5000 

0.3000 

1.0000 

0.5000 

0.3000 

1.0000 

0.5000 

0.3500 

1.0000 

-0.5000 

0.3500 

1.0000 

-0.5000 

0.4000 

1.0000 

0.5000 

0.4000 

1.0000 

0.5000 

0.4500 

1.0000 

-0.5000 

0.4500 

1.0000 

-0.5000 

0.5000 

1.0000 

0.5000 

0.5000 

1.0000 

0.5000 

0.5500 

1.0000 

-0.5000 

0.5500 

1.0000 

-0.5000 

0.6000 

1.0000 

1.0000 

0.5000 

0.6500 

1.0000 

-0.5000 

0.6500 

1.0000 

-0.5000 

0.7000 

1.0000 

0.5000 

0.7000 

1.0000 

0.5000 

0.7500 

1.0000 

1.0000 

-0.5000 

0.8000 

1.0000 

0.5000 

0.8000 

1.0000 

0.5000 

0.8500 

1.0000 

1.0000 

-0.5000 

0.9000 

1.0000 

0.5000 

0.9000 

1.0000 

0.5000 

0.9500 

1.0000 

-0.5000 

0.9500 

1.0000 

-0.5000 

1.0000 

1.0000 

0.5000 

1.0000 

1.0000 

0.5000 

1.0500 

1.0000 

-0.5000 

1.0500 

1.0000 

-0.5000 

1.1000 

1.0000 

0.5000 

1.1000 

1.0000 

0.5000 

1.1500 

1.0000 

-0.5000 

1.1500 

1.0000 

-0.5000 

1.2000 

1.0000 

0.5000 

1.2000 

1.0000 

0.5000 

1.2500 

1.0000 

-0.5000 

1.2500 

1.0000 

-0.5000 

1.3000 

1.0000 

0.5000 

1.3000 

1.0000 

0.5000 

1.3500 

1.0000 

-0.5000 

1.3500 

1.0000 

-0.5000 

1.4000 

1.0000 

0.5000 

1.4000 

1.0000 

0.5000 

1.4500 

1.0000 

-0.5000 

1.4500 

1.0000 

-0.5000 

1.5000 

1.0000 

0.5000 

1.5000 

1.0000 

0.5000 
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Table  4  Table  of  true  values  of  ai,  02,  xq,  xi,  and  X2  with  corresponding  estimates 


ai 

a2 

xo 

Xl 

X2 

A 

ai 

A 

a2 

A 

Xo 

A 

Xl 

X2 

0.2000 

0.3000 

1.0000 

0.5000 

0.3000 

0.2000 

0.2010 

1.0000 

-29.0710 

29.8710 

0.2000 

0.5000 

0.3000 

0.2001 

0.4001 

1.0000 

0.5005 

0.2995 

0.3000 

0.4000 

1.0000 

0.5000 

0.3000 

0.3017 

0.4031 

1.0000 

0.5175 

0.2825 

0.4000 

0.5000 

1.0000 

0.5000 

0.3000 

0.3755 

0.4823 

1.0000 

0.3354 

0.4646 

0.5000 

0.6000 

1.0000 

0.5000 

0.3000 

0.5000 

0.5011 

1.0000 

-27.1234 

27.9234 

0.6000 

0.7000 

1.0000 

0.5000 

0.3000 

0.6500 

0.6510 

1.0000 

10.3104 

-9.5104 

0.6000 

0.8000 

1.0000 

0.5000 

0.3000 

0.6083 

0.8054 

1.0000 

0.5292 

0.2708 

1.0000 

0.5000 

0.3000 

0.6697 

0.7807 

1.0000 

0.3113 

0.4886 

0.7000 

0.9000 

1.0000 

0.5000 

0.3000 

0.6621 

0.8373 

1.0000 

0.2846 

0.5154 

1.0000 

0.5000 

0.3000 

0.8000 

0.8011 

1.0000 

-27.1771 

27.9771 

1.0000 

0.7932 

0.9904 

1.0000 

0.4681 

0.3319 

0.9000 

1.0000 

1.0000 

0.5000 

0.3000 

0.9499 

0.9509 

1.0000 

10.1132 

-9.3132 

1.0000 

1.0031 

11062 

1.0000 

0.5332 

0.2668 

1.0000 

0.5000 

0.3000 

1.0500 

1.1773 

1.0000 

0.2500 

0.5500 

1.2000 

0.5000 

0.3000 

1.2511 

1.2521 

1.0000 

11.5485 

-10.7485 

1.0000 

0.5000 

0.3000 

1.2000 

1.3714 

1.0000 

0.1659 

0.6325 

1.4000 

1.5000 

1.0000 

0.5000 

0.3000 

1.3862 

1.4867 

1.0000 

0.3917 

0.4083 

IV.  Analysis 


4.1  Overview 

The  theory  of  operations  of  the  eMRDLL  and  the  NCDLL  provided  in  Chapters  2 
and  3  describe  the  expected  performance  characteristics  of  these  designs.  However,  ideal 
conditions  were  assumed  in  the  derivations  of  the  principle  concepts  of  these  designs.  As 
such,  analysis  of  these  designs  in  a  realistic  environment  is  necessary  in  order  to  test  their 
effectiveness.  Computer  modeled  simulations  provide  this  realistic  environment. 

This  Chapter  analyzes  the  operation  of  the  eMRDLL  and  the  NCDLL  in  a  simulated 
environment.  The  eMRDLL  is  used  as  the  platform  to  test  the  eflFectiveness  of  ML  estima¬ 
tion  in  a  GPS  environment.  The  NCDLL,  with  correlator  spaicings  of  A  =  1.0  and  A  =  0.1, 
serves  as  the  baseline  design  for  which  the  performance  of  MRDLL  can  be  compared. 

Simiilations  were  conducted  using  a  graphiced  environment  provided  by  the  Simulink 
Toolbox,  version  1.3  which  requires  the  use  of  Matlab,  version  4.2c  to  provide  computa¬ 
tional  support.  In  addition,  functions  from  the  Signal  Processing  Toolbox,  the  Optimiza¬ 
tion  Toolbox,  and  the  Communications  Toolbox  were  used.  All  software  packages  used  in 
these  simulations  were  provided  by  The  Math  Works,  Ihc.  of  Natick,  Massachusetts.  A  de¬ 
scription  of  the  models  developed  is  provided  in  Appendix  A.  Simulations  were  performed 
on  the  Sun  Workstations  provided  at  the  Air  Force  Institute  of  Technology  (AFIT). 

4.2  Simulations  Performed 

Analysis  of  five  simulation  scenarios  are  presented  in  order  to  characterize  the  per¬ 
formance  characteristics  of  eMRDLL. 

•  Simulation  #  1:  Performance  cinalysis  is  conducted  on  the  MCEU  (without 
AWGN)  for  various  multipath  delays  a  €  (0,1.5].  Perfect  synchronization  is  assumed  so 
that  the  MCEU  is  isolated.  Analysis  of  the  MCEU  is  provided  with  and  without  the 
baseband  mixing  process  implemented. 

•  Simulation  #  2:  Performance  analysis  is  conducted  on  the  MCEU  for  vEirious 
multipath  delays  a  G  (0, 1.5]  for  a  range  of  SNRG  [—35, 30]  dB.  Perfect  synchronization  is 
assumed  so  that  the  MCEU  is  isolated. 
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•  Simulation  #  3:  Closed  loop  code  phase  tracking  analysis  is  performed  (with¬ 
out  AWGN)  for  the  eMRDLL  versus  the  NCDLL  (A  =  1.0  and  A  =  0.1)  for  various 
multipath  delays  a  G  (0,1.5]. 

•  Simulation  #  4:  Closed  loop  code  phase  tracking  analysis  is  performed  for 
the  eMRDLL  for  various  multipath  delays  a  G  (0,1.5]  for  a  range  of  SNRg  [—35,30]  dB. 

•  Simulation  ^  5:  Performance  analysis  is  conducted  on  the  MCEU  for  time 
varying  multipath  delays  with  and  without  AWGN.  Perfect  synchronization  is  assumed  so 
that  the  MCEU  is  isolated. 

4.3  Simulation  Parameters 

In  order  to  provide  an  accurate  analysis  of  the  eMRDLL  and  the  NCDLL,  signal 
parameters  were  chosen  so  as  to  best  emulate  the  functionality  of  a  receiver  in  a  GPS 
environment.  Software  limitations,  however,  existed  such  that  meiny  parameters  foimd  in 
GPS  applications  had  to  be  scaled.  The  following  sections  outline  the  parameters  used  for 
simulations  conducted  as  well  as  the  rationale  for  choosing  these  parameters. 

In  Choosing  the  simulation  parameters  properly,  the  following  were  considered  as 
constraints; 

•  In  a  simulation  environment,  such  as  Simulink,  the  sampling  rate,  /,,  must 
be  finite.  Increasing  the  sampling  rate  increases  simulation  run-times,  so  there  exists  a 
trade-off  between  computational  burden  and  aliasing  considerations. 

•  The  code  rate,  1/Tc  must  be  chosen  high  enough  to  ensure  proper  operation 
of  receiver  components  such  as  the  VCC;  for  example,  a  negative  input  into  the  VCC  with 
a  code  rate  that  is  very  low  may  drive  the  frequency  of  the  VCC  to  zero  at  which  point  it 
cannot  be  lowered  further  (i.e,  cannot  produce  negative  frequency). 

•  Filters  with  very  narrow  bandwidths,  implemented  in  Simulink,  require  a 
carrier  frequency  to  code  rate  ratio,  fo  :  1/Ic,  not  to  exceed  10:1  in  order  to  ensure  proper 
performance.  Paxeuneters  chosen  which  exceed  this  ratio  lead  to  amplification  of  the  signal 
within  the  filter’s  passband. 
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4.3.1  Spreading  Code  Chip  Rate,  1/Tc.  The  determination  of  the  chip  rate,  1/Tc, 
was  a  fundamental  consideration  in  designing  the  models  of  the  eMRDLL  and  the  NCDLL. 
A  chip  rate  of  1/Tc  =  100  chips/ sec  (or  Hz)  Wcis  chosen  to  ensure  proper  operation  of 
receiver  components  such  as  the  VCC.  The  clock  rate  of  the  VCC  had  to  be  considered 
such  that  it  would  never  reach  zero  (the  lower  limit  of  the  VCC  dyneunic  range)  based  on 
the  error  signal  coming  from  the  loop  filter,  F(s).  Depending  on  the  gain  of  the  system,  a 
chip  rate  less  than  100  chips/ sec  could  result  in  the  VCC  being  driven  to  zero.  A  higher 
chip  rate  also  allowed  for  more  fiexibility  in  choosing  the  loop  bandwidth,  Bl,  for  the 
NCDLL  (i.e.,  smaller  loop  bandwidths  could  be  chosen  for  higher  chip  rates). 

4.3.2  Carrier  Frequency,  fo  .  In  order  to  ensure  that  the  filters  performed 
properly,  a  carrier  frequency  of  /©  =  1000  Hz  was  chosen  so  as  to  maintain  a  carrier  to 
code  chip  rate  ratio  of  10:1.  This  also  helped  keep  the  sampling  frequency,  lower  which 
was  necessary  to  maintain  reasonable  simulation  run  times. 

There  was,  however,  a  trade-off  in  choosing  a  carrier  frequency  which  was  relatively 
close  to  the  code  rate.  As  will  be  shown  in  the  following  sections,  higher  frequency  har¬ 
monics,  which  resulted  from  the  pre-correlation  mixing  process,  were  introduced  into  the 
baseband  signal  entering  the  eMRDLL. 

4.3.3  Sampling  Frequency,  fs.  The  sampling  frequency,  Wcis  chosen  rela¬ 
tive  to  the  carrier  frequency,  fo,  and  the  code  chip  rate,  l/Tc-  A  chip  resolution  of  100 
samples/chip  was  desired  in  order  to  be  able  to  consider  a  broad  range  of  multipath  de¬ 
lays.  Resolution  less  than  .OlTc  would  have  resulted  in  a  more  limited  range  of  possible 
multipath  delays.  Also,  at  least  10  samples  per  carrier  cycle  were  desired  for  Nyquist 
rate  considerations  regarding  representation  of  the  received  multipath  signed.  In  order  to 
maintain  reasonable  simulation  nm  times  and  meet  the  sampling  requirements  described 
above,  fs  was  chosen  as  10,000  Hz.  This  resulted  in  a  simulation  integration  step  size  of 
l/fs  =  0.0001  sec. 
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C/A  (DS/SS)  Code  Length,  N.  GPS  satellites  transmit  a  Gold  Code^ 
with  meiximum  length  pseudo-noise  (PN)  characteristics.  As  stated  in  Chapter  1,  the  C/A 
spreading  code  has  N  =  1023  chips  per  code  cycle.  Preliminary  simulations  indicated 
that  the  response  time  of  using  an  JV  =  1023  chip  sequence  was  too  lengthy.  As  such,  a 
maximum  length  PN  sequence  of  iV^  =  2”*  —  1  =  63  chips,  with  code  shift  register  order 
of  m  =  6  was  used,  as  described  in  Appendix  A.  Choosing  AT  =  63  for  the  simulations 
still  permits  the  use  of  the  large  N  code  correlation  approximation  of  Equation  8  with  the 
added  benefit  of  increasing  simulation  performance. 

4.3.5  BPF  and  LPF  Bandwidths,  B,  As  stated  in  Chapter  2,  the  bandwidths 
of  the  filters  in  the  ‘early’-'late’  gate  branches  should  be  chosen  on  the  order  of  the  data 
rate.  In  GPS,  this  corresponds  to  50  if z  which  is  a  fraction  of  the  spreading  code  chip 
rate  of  1.023  MHz.  Trying  to  set  the  filter  bandwidths  of  the  models  used  in  Simulink 
with  the  same  proportion  of  that  in  GPS  would  have  resulted  in  imrealizable  constraints 
placed  on  the  filters,  van  Nee  states  that,  for  GPS,  these  pre-detection  bandwidths  are 
typically  in  the  range  of  1  kHz  down  to  the  two-sided  null-to-null  data  bandwidth  of  100 
Hz  [19].  The  integration  period  determined  from  the  pre-detection  bandwidth  of  1  kHz 
would  be  equivcdent  to  1/NTc.  This  is  more  reeisonable  for  a  simulation  environment  than 
a  bandwidth  chosen  on  the  order  of  the  data  rate. 

The  code  correlation  function,  Rc{t),  of  Equation  7,  can  be  implemented  with  a 
code  multiplier  followed  by  either  a  BPF  with  bandwidth  B  =  2/NTc  Hz  or  LPF  with 
bandwidth  B  =  1/NTc,  where  either  filter  acts  as  an  integrator.  As  such,  the  LPFs  used  in 
eMRDLL  were  set  at  B  =  IfNTc  =  1.6  Hz.  In  the  NCDLL  models,  the  BPF  bandwidths 
could  not  be  set  at  B  =  2/NTc  =  3.2  Hz  due  to  filter  bandwidth  constraints  encountered 
with  Simulink;^  consequently,  the  bandwidths  of  the  BPFs  had  to  be  set  to  a  slightly 
higher  value  of  B  =  4  Hz. 

^The  gold  code  represents  the  modulo- two  addition  of  two  maximum  length  sequences  generated  by  10 
bit  registers. 

^Filter  bandwidths  that  were  too  narrow  resulted  in  signals  that  were  improperly  filtered. 
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^.5.6  Loop  Filter  Bandwidth,  Bi.  In  most  GPS  receivers,  the  loop  integration 
time  is  roughly  the  reciprocal  of  the  code-tracking  loop  bandwidth.  Bandwidths  of  code 
tracking  loops  can  vary  between  0.02  Hz  and  1  Hz  [19,21].  In  conducting  the  simulations, 
the  loop  filter  bandwidth  could  not  be  scaled  lower,  because  DC  terms  would  have  been 
attenuated  in  the  same  manner  as  higher  firequency  components.  The  first  order  ax;tive 
lead-lag  filter  should  integrate  DC  terms;  therefore,  a  loop  filter  bandwidth  was  chosen 
within  the  range  of  an  actued  GPS  receiver.  A  loop  filter  bandwidth  of  Bi  =  0.1  Hz 
(wn  =  O.OGtt  rad! sec  and  C  =  l/\/2  in  by  Equation  89)  provided  the  best  performance 
for  the  NCDLL  model  implemented  in  Simulink.  The  parameters  a;„  =  O.OGtt  rad/ sec 
and  C  =  l/\/2  were  left  unchanged,  relative  to  the  NCDLL,  in  designing  the  ALC  of  the 
eMRDLL;  however,  the  loop  bandwidth  does  vary  with  A  due  to  the  presence  of  the  ALC. 

4.S.  7  Signal  to  Noise  Ratio,  SNR  .  SNR  was  computed  via  the  ratio  of  the  power 
of  the  transmitted  signal,  Poq,  to  the  noise  power,  N,  which  is  determined  by  the  variance 
of  the  noise,  <t^.  In  using  the  AWGN  noise  block,  provided  by  the  Communications  Toolbox 
for  use  in  Simulink,  it  was  discovered  that  the  sampled  noise  output  had  a  different  variance 
over  a  period  of  100,001  samples  than  what  was  specified  as  the  input  variemce.  This  was 
true  for  all  variance  inputs  used  in  the  simulations.  As  such  the  SNR  levels  specified  did 
not  correspond  with  calculated  SNR  levels.  Table  5  provides  the  input  SNR  versus  the 
calculated  SNR  used  in  the  simulations  discussed  in  this  chapter. 


Table  5  Specified  SNR  versus  calculated  SNR 


Specified  SNR  (dB) 

-35 

-23 

-13 

-6 

0 

Q 

10 

15 

20 

25 

30 

Calculated  SNR  (dB) 

-8 

-4.5 

a 

3.5 

6 

8.5 

11 

13.5 

4-3.8  Simulation  Parameter  Summary  .  The  following  provides  a  list  of  all 
simulation  parameters  used  in  the  following  analysis.  Parameters  described  below  which 
were  not  discussed  in  the  previous  sections  were  based  on  the  designer’s  preference. 

•  Code  chip  rate:  l/Tc  =  100  chips/ sec  (or  100  Hz) 

•  Carrier  frequency:  /©  =  1000  Hz 
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•  Simulation  sampling  frequency:  /,  =  10000  Hz 

•  Code  length:  N  =  6Z 

•  BPF  bandwidth  (NCDLL):  B  =  i  Hz 

•  BPF  type:  2”'^  Order  Butterworth 

•  LPF  bandwidth  (eMRDLL  and  NCDLL):  B  =  1.6  Hz 

•  LPF  type:  5*^  Order  Butterworth 

•  Loop  filter  bandwidth:  0.1  Hz 

•  TYansmitted  signal  power:  P  =  1/2 

•  Signal  amplitudes:  oq  =  1.0  and  ai  =  0.5 

•  Direct-path  propagation  delay:  To  =  0 

•  ‘Early-late’  normalized  correlator  spax:ing  (NCDLL):  A  =  1.0,  A  =  0.1 

•  ‘Early-late’  normalized  correlator  spacing  (eMRDLL):  A  =  0.1 

•  Loop  natural  frequency:  a;»  «  O.OStt  rad/ sec 

•  Loop  damping  factor:  (  =  l/\/2 

•  Simulink  (relative)  SNR  levels:  -35  dB  to  30  dB 

44  Simulation  #  1:  MCEU  Noise  Free  Performance 

The  code  phase  tracking  accuracy  of  the  eMRDLL  is  based  largely  on  the  perfor¬ 
mance  of  the  MCEU.  As  such,  the  MCEU  was  analyzed  to  characterize  its  behavior  and 
performance  characteristics.  The  analysis  in  the  following  section  was  performed  assuming 
that  the  MCTL  was  perfectly  synchronized  at  all  times,  which  isolated  the  performance  of 
the  MCEU.  This  analysis  was  performed  without  AWGN;  a  noise-free  case. 

44’^  MCEU  Performance  (without  baseband  mixing  operation).  As  discussed  in 
section  4.3.2,  the  choice  of  the  carrier  frequency,  /o,  resulted  in  some  harmonic  spillage  into 
the  spectrum  of  the  baseband  signal.  The  objective  of  this  anedysis  is  to  charaxiterize  the 
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performance  of  the  MCEU  for  the  case  where  the  spreading  code  is  already  at  baseband; 
thus,  removing  the  effects  of  the  frequency  conversion  from  IF  to  baseband. 

In  implementing  the  MCEU,  the  rate  at  which  the  MCEU  provides  estimates  to 
the  MCTL  and  ALC  must  be  chosen.  For  a  system  operating  at  real-time,  it  would  be 
undesirable  and  unnecessary  to  estimate  the  multipath  parameters  at  a  sampling  rate 
higher  than  the  code  rate;  a  practical  system  requires  estimates  to  be  provided  at  most 
once  every  code  period,  iV'Tc,  depending  on  the  application  of  the  receiver.^  The  MCEU 
was  designed  to  provide  estimates  to  the  MCEU  once  every  code  period,  NTc  =  0.63  sec. 
Since  filters  are  being  used  as  integrators,  the  absolute  time  at  which  estimates  are  sampled 
is  very  important.  To  demonstrate  the  importance  of  sampling  time,  the  estimates  from 
the  MCEU  were  calculated  every  one-tenth  of  a  cycle,  0.063  sec.  Figures  37  and  38  provide 
plots  of  the  estimates  of  d,  io)  and  xi.  Each  figure  illustrates  the  true  values  of  a,  a?0j  and 
xi  as  well  as  estimates  of  these  parameters  at  sampling  intervals  of  0.63  sec  and  0.063  sec. 
The  choice  of  where  sampling  occurred  was  empirically  determined  based  on  the  response 
of  the  LPFs  implemented. 


a  =  0.5 


Figure  37  MCEU  performance  without  baseband  conversion  mixing  process  (assuming  perfect 
tracking)  where  represents  the  actual  value  of  a,  ‘ — ’  represents  MCEU  estimates 
of  a  sampled  every  one-tenth  of  a  cycle,  and  -’represents  MCEU  estimates  of  a 
sampled  once  every  cycle. 


^Receivers  operating  in  environments  with  slowly  time  varying  multipath  would  not  need  to  update 
signal  estimates  as  frequently  as  receivers  operating  in  a  dynamic  multipath  environment. 
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a  =  0.5  a  =  0.5 


Figure  38  MCEU  performance  without  baseband  conversion  mixing  process  (assuming  perfect 
tracking)  where  represents  the  actual  value  of  Xo  and  xi,  ‘ — ’  represents  MCEU 
estimates  of  xo  and  xi  sampled  every  one-tenth  of  a  cycle,  and  -’represents  MCEU 
estimates  of  xq  and  xi  sampled  once  every  cycle. 

From  the  figures,  the  seunpling  interval  of  0.63  sec  is  evident  as  the  period  of  the 
steady-state  response  of  the  estimates  provided  by  the  MCEU  for  a,  xq,  and  ®i;  these 
oscillations  are  caused  by  the  LPFs  used  in  the  bank  of  correlators.  The  code  autocorre¬ 
lation,  E(t),  function  provided  in  Equation  7,  indicates  that  R(t)  is  time  averaged  over 
a  single  period.  One  can  consider  the  operation  of  R{t)  over  time  as  a  moving  window 
integrator  that  calculates  the  area  of  the  received  spreeiding  code  coinciding  with  a  locally 
generated  code  replica,  delayed  by  l3kTc,  for  each  corresponding  correlator  ^lrm.  Assuming 
perfect  synchronization,  R{t)  would  remain  constant  for  an  ideal  integrator;  however,  the 
LPFs  used  in  the  correlator  arms  are  not  ideal.  The  outputs  of  the  correlator  arms  vary 
with  time  as  shown  in  the  figmes.  This  results  in  the  MCEU  producing  estimates  that  also 
vary  with  time.  As  it  was  noted,  the  responses  of  the  filters  and  thus  the  estimates  of  the 
MCEU  are  periodic  with  period  equal  to  the  length  of  the  code  period,  NTc]  consequently, 
the  figures  illustrate  that  the  MCEU  estimates  indeed  have  a  period  of  NTc  =  0.63  sec. 

Based  on  the  results  provided  above,  one  can  time  estimator  performance  by  adjust¬ 
ing  the  epoch  at  which  the  MCEU  periodic  samples  are  collected;  this  can  be  established 
a  priori  based  on  the  response  characteristics  of  the  particular  LPF  implemented. 
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Figures  37  and  38  indicate  that  the  the  saimpling  time  implemented  for  the  sampling 
rate  of  1/.63,  was  chosen  very  well  since  the  estimates,  a  and  ®0)  are  very  near  to  the  true 
values,  a  and  xq.  The  estimate  xi,  however,  is  not  as  accurate  as  the  estimates  of  the 
other  two  parameters. 

The  sensitivity  of  ®i  can  best  be  described  by  considering  that  d  ^  a  causes  ^  4>e 
and  the  baseband  amplitudes  of  ajQ  and  aji  deviate  as  presented  in  Figure  39;  this  figure 
provides  the  expected  values  of  so  and  xi  for  signed  amplitudes  of  oq  =  1.0  and  oi  =  .5 
as  a  is  varied  through  (0,1.5].  It  is  interesting  to  note  that  the  range  of  amplitude  of  the 
reflected  component  is  much  larger  than  that  of  so  and  that  the  envelope  amplitude  of  so 
exponentially  decreases  to  zero  as  a  increases  from  zero.  Figure  39  thus  illustrates  how 
the  estimates  sq  and  si  can  be  expected  to  oscillate  with  oscillating  values  of  d. 


a 


Figure  39  Baseband  values  of  sq  and  sj  for  P  =  1/2  and  /©Tc  =  10. 

4.4-^  MCEU  Performance  (with  baseband  mixing  operation).  Now  that  the 
baseline  performance  of  the  MCEU  has  been  char«icterized,  the  effects  of  the  baseband 
conversion  process  on  the  MCEU  can  be  analyzed.  Figures  40  through  45  illustrate  the 
performance  of  the  MCEU  for  three  cases,  where  a  =  0.1, 0.5,  and  1.3.  As  with  the  previ¬ 
ous  section,  the  estimates  were  calculated  every  0.063  sec  as  well  as  estimates  at  sampling 
intervals  of  0.63  sec.  Figures  axe  provided  for  d  and  »o»  »1  for  each  of  the  multipath  cases 
examined  in  this  section. 
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4.4.2. 1  MCEU  Performance  for  a  =  0.1.  The  estimates  of  a,  xq,  and 
XI,  provided  in  Figures  40  and  41,  illustrate  the  sensitivity  caused  by  the  LPF  responses 
at  small  delays  of  a,  such  as  a  =  0.1.  The  estimates  provided  by  the  MCEU  cover  the 
fuU  range  of  possible  values  of  a.  The  estimates  are  also  periodic  in  nature  with  period 
NTc  as  previously  described.  Although  the  samples  taken  every  0.63  seconds  provide  an 
estimate  a  =  0.15,  most  0.063  sec  estimates  are  far  away  from  the  true  a  =  0.1.  The 
MCEU  estimates  »o  and  ®i  also  vary  greatly  based  on  varying  a  estimates.  As  presented 
in  Figure  39,  both  jcq  and  x\  are  very  sensitive  for  variations  of  a  near  the  direct-path 
signal.  Figure  41b  also  shows  that  as  or  varies,  xi  may  become  negative. 


Figure  40  MCEU  performance  with  baseband  conversion  mixing  process  (assuming  perfect 
tracking)  where  represents  the  actual  value  of  a,  ‘ — ’  represents  MCEU  esti¬ 
mates  of  a  sampled  every  one-tenth  of  a  cycle,  and  -’represents  MCEU  estimates 
of  a  sampled  once  every  cycle. 

4.4-2. 2  MCEU  Performance  for  a  =  0.5.  Figures  42  and  43  illustrate 
estimates  of  a,  xq,  and  for  a  =  0.5.  The  responses  of  the  MCEU,  for  the  diflferent 
sampling  rates,  are  identical  to  those  described  in  the  section  where  baseband  mixing 
process  was  neglected  (with  the  exception  of  the  amplitudes  of  ®o  and  xi).  Notice  that 
the  amplitudes  of  xq  and  xi  are  less  than  the  true  amplitudes  of  xq  and  xi.  This  is  a 
direct  result  of  baseband  mixing  process.  The  decrease  in  estimated  values  are  attributed 
to  the  low  ratio  of  fo  :  1/Tc  in  which  harmonics  of  the  higher  frequency  terms  are  not 
low  pass  filtered.  These  harmonics  decrease  the  level  of  energy  present  in  the  correlator 
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Figure  41  MCEU  performance  with  baseband  conversion  mixing  process  (assuming  perfect 
tracking)  where  represents  the  actual  value  of  xq  and  xi,  represents  MCEU 
estimates  of  xq  and  xi  sampled  every  one-tenth  of  a  cycle,  and  -’represents  MCEU 
estimates  of  xo  and  xi  sampled  once  every  cycle. 


samples  resulting  in  estimates  of  xq  <  ®o  and  xi  <  xi.  Since  the  correlator  measurements 
of  each  arm,  maintain  a  fixed  ratio  with  respect  to  each  other,  &  is  not  degraded.  Smaller 
estimates  of  xo  and  xi,  relative  to  true  parameters,  occurred  for  all  cases  of  a  tested;  the 
results  of  these  test  cases  are  summarized  in  Section  4.4.3. 

It  should  be  noted  that  underestimating  the  values  of  xq  and  xi  should  not  present  a 
significant  problem  in  the  code  tracking  performance  of  the  MCTL  since  the  same  baseband 
conversion  process  applies  to  the  MCTL.  As  such,  the  relative  amplitudes  of  the  signals 
processed  in  the  MCTL  will  be  proportional  to  the  amplitudes  of  the  signals  processed  in 
the  MCEU. 


4.4-^’3  MCEU  Performance  for  a  =  1.3.  Figures  44  and  45  provide  very 
similar  results  to  those  previously  described.  A  significant  difference  is  the  estimated  value 
d,  illustrated  in  Figure  44,  which  has  larger  estimation  error  relative  to  the  estimates 
d  €  [.2,1.0].  This  increased  estimation  error  is  present  for  d  6  (1.0, 1.5].  Note  that  at 
the  sample  rate  of  1/.63  Hz,  estimates  are  collected  from  the  lower  bound  of  the  response 
curve  envelope.  The  true  value  of  a  is  located  at  the  center  of  the  MCEU  response  curve 
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a  =  .5 


Figure  42  MCEU  performance  with  baseband  conversion  mixing  process  (assuming  perfect 
tracking)  where  represents  the  actual  value  of  a,  ’  represents  MCEU  esti¬ 
mates  of  a  sampled  every  one-tenth  of  a  cycle,  and  -’represents  MCEU  estimates 
of  a  sampled  once  every  cycle. 

envelope.  For  this  case,  the  epochs  at  which  samples  were  teiken  led  to  a  less  accurate  d  , 
relative  to  the  a  =  0.5  case. 

4,4.3  Summary  of  MCEU  Performance  .  The  MCEU  was  tested  for  values  of  a 
at  increments  of  0.1  through  the  range  a  €  [0.1, 1.5],  with  and  without  the  beiseband  mixing 
operation  (BBMO)  implemented.  Based  on  the  results  summarized  in  Table  6,  Figure  46 
provides  the  estimates  of  a  and  Figure  47  the  estimates  of  xq  and  afi,  each  versus  the  actual 
values.  Figure  46  illustrates  the  accuracy  potential  of  the  MCEU.  All  estimates  with  the 
exception  the  estimate  of  a  =  0.1  are  within  0.03  chips  of  the  true  Veilue  of  a.  The  figure 
also  illustrates,  as  previously  discussed,  that  the  baseband  conversion  process  does  not 
affect  the  estimates  of  a.  Figure  47,  however,  indicates  that  the  baseband  conversion 
process  does  affect  the  estimates  of  xq  and  xi.  Furthermore,  with  the  exception  of  the 
estimate  for  a  =  0.1,  the  difference  between  the  estimates  of  »o  and  ®l  with  baseband 
conversion  versus  no  baseband  conversion  is  nearly  constant  throughout  the  range  of  a 
tested.  In  summary,  d  is  determined  by  the  relative  differences  between  the  values  of  iJjt; 
whereas,  xq  and  ®i  are  directly  affected  by  the  amplitudes/energy  levels  at  the  correlator 
inputs. 
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a  =  .5  a  =  .5 


Figure  43  MCEU  performance  with  baseband  conversion  mixing  process  (assuming  perfect 
tracking)  where  represents  the  actual  value  of  a;o  and  Xi,  ‘ ’  represents  MCEU 
estimates  of  Xq  and  xi  sampled  every  one-tenth  of  a  cycle,  and  -’represents  MCEU 
estimates  of  xo  and  xi  sampled  once  every  cycle. 


a  =  1.3 


Figure  44  MCEU  performance  with  baseband  conversion  mixing  process  (assuming  perfect 
tracking)  where  represents  the  actual  value  of  a,  ‘ — ’  represents  MCEU  esti¬ 
mates  of  a  sampled  every  one-tenth  of  a  cycle,  and  -’represents  MCEU  estimates 
of  a  sampled  once  every  cycle. 
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a  =  1.3  a  =  1.3 


Figure  45  MCEU  performance  with  baseband  conversion  mixing  process  (assuming  perfect 
tracking)  where  represents  the  actual  value  of  xo  and  a;i,  ‘ ’  represents  MCEU 
estimates  of  xq  and  xi  sampled  every  one-tenth  of  a  cycle,  and  -’represents  MCEU 
estimates  of  xq  and  Xi  sampled  once  every  cycle. 


Table  6  Summary  of  MCEU  performance  where  a,  xo  =  1.0,  and  xi  =  0.5  are  the  true  param¬ 
eters. 


with  BBMO 

without  BBMO 

a 

a 

XQ 

XI 

a 

Xo 

^1 

0.1000 

0.1500 

1.0883 

0.3161 

0.1500 

1.1633 

0.3380 

0.2000 

0.1974 

0.9332 

0.4701 

0.1974 

0.9976 

0.5026 

0.3000 

0.3000 

0.9346 

0.4708 

0.3000 

0.9990 

0.5033 

0.4000 

0.4013 

0.9362 

0.4699 

0.4013 

1.0007 

0.5024 

0.5000 

0.5000 

0.9328 

0.4762 

0.5000 

0.9971 

0.5091 

0.6000 

0.5926 

0.9300 

0.4766 

0.5926 

0.9941 

0.5096 

0.7000 

0.6997 

0.9360 

0.4709 

0.6997 

1.0005 

0.5034 

0.8000 

0.7995 

0.9361 

0.4704 

0.7995 

1.0007 

0.5029 

0.9000 

0.8994 

0.9361 

0.4704 

0.8994 

1.0007 

0.5029 

1.0000 

0.9912 

0.9323 

0.4659 

0.9912 

0.9967 

0.4981 

1.1000 

1.0758 

0.9296 

0.4507 

1.0758 

0.9937 

0.4819 

1.2000 

1.1881 

0.9228 

0.4476 

1.1881 

0.9865 

0.4786 

1.3000 

1.2785 

0.9160 

0.4286 

1.2785 

0.9792 

0.4583 

1.4000 

1.3736 

0.9092 

0.4078 

1.3736 

0.9720 

0.4361 

1.5000 

1.4892 

0.9024 

0.4061 

1.4891 

0.9647 

0.4342 
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Figure  46  MCEU  estimates  of  a  versus  actual  values  of  a  (MCEU  sampling  period:  0.63  sec). 

Note:  ‘o’  represents  MCEU  estimates  with  mixing  operation;  ‘x’  represents  MCEU 
estimates  without  mixing  operation.  (MCEU  sampling  period:  0.63  sec). 


0.2  0.4  0.6  0.8  1  1.2  1.4 


a 


Figure  47  Estimated  values  of  xq  and  xi.  Note:  ‘o’  represents  MCEU  estimates  with  mixing  op¬ 
eration;  ‘x’  represents  MCEU  estimates  without  mixing  operation.  (MCEU  sampling 
period:  0.63  sec). 
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4.5  Simulation  #  2;  MCEU  Performance  in  AWGN 

Since  the  performance  baseline  of  the  MCEU  has  been  characterized  imder  ideal  con¬ 
ditions,  the  MCEU  performance  will  now  be  examined  in  the  presence  of  AWGN.  In  partic¬ 
ular,  it  is  desired  to  determine  the  SNR  threshold  at  which  the  MCEU  performance  begins 
to  significantly  degrade.  SNR  levels  were  varied  through  the  range  of  SNRG  [—35,30]  dB 
as  specified  in  Table  5.  Once  again,  the  MCTL  is  assumed  perfectly  synchronized  with  the 
received  signal.  MCEU  samples  are  output  every  0.63  sec  over  a  period  of  10  sec,  providing 
15  estimates  for  analysis. 

Figure  48  illustrates  the  meein  of  the  estimates  of  a,  ®o>  and  ®i  for  delays  of  a  € 
[0.1, 1.5]  in  increments  of  0.1  (NOTE:  These  figures  are  only  presented  to  show  trends, 
not  to  provide  detail).  Each  line  of  Figure  48a  represents  the  mean  of  the  estimates  of 
an  a  G  [0.1, 1,5].  The  figure  indicates  that  the  mean  of  a  and  io>  ®l  are  approximately 
constant  for  SNRG  [—23, 30]  dB  with  the  exception  of  the  estimates  of  delays  a  =  0.1  and 
a  =  0.2.  The  mean  of  the  estimates  at  these  delays  show  degradation  at  SNR<  0  dB  and 
SNR<  — 13  dB,  respectively. 


1.6| - 1 - 1 - 1 - 1 - 1 - 1 - 1  I - 1 - 1 - r 


"  -’§NR(<!b)  -’gNR(^B)  “  = 

a)  Mean  of  a.  b)  Mean  of  xq  and  xi. 

Figure  48  Mean  of  estimates  a,  xq,  and  xi  provided  by  MCEU  for  delays  a  G  [0.1, 1.5]. 


Figure  49  illustrates  the  variance  of  the  estimates  of  a,  xq,  and  xi  for  delays  of 
a  G  [0.1, 1.5]  in  increments  of  0.1  (NOTE:  These  figures  are  only  presented  to  show  trends, 
not  to  provide  detail).  Variance  degradation  is  significant  for  SNR<  —23  dB  for  all  a. 
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Notice  that  the  change  in  the  slopes  of  the  variance  of  the  estimates  for  SNR<  —23  dB 
differ  considerably  through  the  range  of  delays.  The  veiriance  degradation  for  a  =  0.1  and 
a  =  0.2  begins  at  0  dB  and  -13  dB,  respectively. 


Figure  49  Variance  of  estimates  d,  »o>  and  ®l  provided  by  MCEU  for  delays  a  G 
[0.1, 1.5]. 


In  order  to  more  closely  examine  the  noise  threshold  effect,  five  delays  (a  =  0.3, 
a  =  0.5,  a  =  0.7,  a  =  1.0,  and  a  =  1.3)  were  examined  in  detail  through  the  full  range  of 
SNRs.  These  delays  represent  a  uniform  sampling  of  the  15  delays  previously  considered. 
For  these  delays,  the  root  mean  square  error  (rmse)  of  15  samples  taken  over  a  10  second 
window  is  used  as  the  measure  of  performance.  The  rmse  refiects  both  the  mean  error  as 
well  as  the  variance  of  the  estimates. 

Figures  50  and  51  illustrate  the  rmse  of  d,  ®o>  ®nd  «i  for  the  five  delays  examined. 
As  with  the  mean  and  the  variance  previously  discussed,  significant  degradation  in  per¬ 
formance  occurs  for  SNR<  —23  dB.  In  addition,  note  that  little  or  no  noise  effects  are 
observed  when  SNR>  0  dB.  The  estimates  d  for  a  =  1.3  are  most  affected  by  increased 
SNR  levels;  however,  this  is  also  true  for  the  noise-free  case  presented  in  Simulation  1. 

Figures  52  and  53  illustrate  the  rmse  of  estimates  for  all  delays  a  G  [0.1, 1,5]  in¬ 
cremented  by  0.1  at  SNR  levels  of  -23  dB,  0  dB,  and  10  dB.  These  figures  illustrate  that 
there  is  a  much  larger  MCEU  performance  degradation  between  0  dB  and  -23  dB  than 
between  10  dB  and  0  dB.  This  is  consistent  with  the  previously  determined  noise  threshold 
of  SNR=-23  dB. 
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Figure  50  rmse  of  a  provided  by  the  MCEU  for  a  =  0.3  (’ — ’),  a  =  0.5  (’-*’)>  =  0.7  (’- 

a  =  1.0  and  a  =  1.3 

As  previously  indicated,  MCEU  performance  is  the  most  noise  sensitive  for  the  small 
reflected  signal  delays  of  a  =  0.1  and  a  =  0.2;  a  for  delay  a  =  1.3  shows  the  most 
sensitivity  to  increased  noise  levels,  d  of  a  =  0.5  is  least  affected  by  increased  noise  levels. 

Other  than  the  estimates  of  xq  and  ®i  for  delays  a  =  0.1  and  a  =  0.2,  ®o  ®i 
for  delay  a  =  0.3  show  the  largest  rmse  for  estimates  of  xq  and  ®i.  Degradation  in  the 
estimates  ®o  and  ®i  is  uniform  for  a  decrease  in  SNR  from  10  dB  to  0  dB‘,  however,  the 
increase  in  rmse  due  to  a  decrease  in  SNR  from  0  dB  to  -23  dB  is  somewhat  random  for 
delays  a  6  [0.3, 1.5]. 

4.5.1  Summary  of  Results.  When  the  MCTL  is  perfectly  tracking  the  direct  path 
code  phase,  the  MCEU  exhibits  a  rmse  noise  threshold  of  -23  dB,  for  reflected  signal  delays 
of  a  G  [0.3, 1.5];  from  Table  5,  -23  dB  corresponds  to  a  computed  SNR  of  -13  dB. 

A  comparison  of  the  results  from  this  analysis  to  GPS  applications  C2in  be  made  by 
normalizing  a  typical  GPS  SNR.  Spilker  demonstrates  that  a  typical  GPS  noise  density 
level  of  No  =  205.2  dBW/Hz  exists  for  an  equivalent  noise  temperature  of  a  typical  low 
noise  amplifier  configuration  [8].  Given  that  the  bandwidth  of  the  main  lobe  of  the  C/A 
spreading  code  is  5  =  2.046  MHz,  an  approximation  of  the  noise  power  at  the  input  of 
a  receiver  antenna  can  be  expressed  as  N  =  NqB  =  —141.89  dBW.  GPS  requirements 
specify  a  minimum  signal  strength  of  —160  dBW  for  C/A  code  on  the  LI  carrier  which 
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Figure  51  rinse  of  xo  and  xi  provided  by  the  MCEU  for  a  =  0.3  (’ — ’),  a  =  0.5  (’-*’))  Q:  =  0.7 
(’-  a  =  1.0  and  a  =  1.3 


would  provide  an  SNRw  —18dB.  Prom  table  5,  a  specified  SNR  of  -35  dB  resulted  in  a 
calculated  SNR  of  -19  dB.  At  this  SNR,  results  from  Simulation  #  2  indicate  significant 
degradation.  Improvements  in  performance  of  the  eMRDLL  would  most  likely  occur  if  the 
bandwidth  of  the  LPF  could  reduce  to  a  normalized  value  proportioned  value  to  the  data 
rate  «  1/NTc.  In  addition,  using  a  code  length  N=1023  may  provide  improvement  by 
increasing  processing  gain. 
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4.6  Simulation  #  3:  eMRDLL  versus  NCDLL  (Noise-Pree  Comparison) 

The  objective  of  this  section  is  to  determine  the  noise-free  performance  of  the  eM¬ 
RDLL  in  comparison  with  the  NCDLL  (normalized  correlator  spacings  of  A  =  1.0  and 
A  =  0.1).  The  normalized  code  phase  error  (To(t)  -  To{t))/Tc  of  the  code  generator  with 
respect  to  the  received  signal  is  used  as  the  measure  of  performance  of  each  design.  Simu¬ 
lations  were  conducted  without  AWGN  in  order  to  characterized  the  performance  of  each 
design  under  ideal  conditions. 

In  order  to  isolate  the  performance  of  the  MCTL,  initial  open  loop  tests  were  con¬ 
ducted  assuming  perfect  synchronization  by  the  VCC  and  perfect  estimates  from  the 
MCEU.  Tests  were  conducted  with  and  without  the  baseband  mixing  operation  imple¬ 
mented.  Results  indicate  that  a  DC  bias  exists  at  the  discriminator  output,  e(5,t),  as 
shown  in  Figure  54,  for  both  cases.  Although,  the  amplitudes  of  each  of  the  errors  signals 
in  the  figure  are  different,  the  bias  present  is  approximately  equal  between  both  discrimi¬ 
nator  outputs.  Therefore,  the  bias  present  is  likely  attributed  to  a  Simulink  computational 
anomaly  in  the  evaluation  of  the  correlation  of  the  received  direct  path  and  refiected  signals 
with  the  locally  generated  code. 


This  bias  has  an  adverse  effect  on  the  closed  loop  tracking  performance  of  eMRDLL. 
The  ALC  operates  similarly  to  a  basic  first  order  loop  filter.  Any  DC  component  is 
integrated  by  the  ALC,  which  drives  the  VCC  away  from  the  desired  operating  point. 
Because  the  bias  is  relatively  small,  a  degradation  in  tracking  performance  does  not  occur 
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for  the  first  several  cycles  in  most  of  the  test  cases.^  The  bias  has  the  least  effect  when 
the  refiected  signal  delay  is  a  =  0.3,  as  shown  in  Figure  55a,  in  which  the  bias  does  not 
impact  synchronization  until  t  23  sec  («  36.5  code  cycles).  The  bias  heis  the  greatest 
effect  when  the  reflected  signal  delay  is  a  =  0.5,  as  shown  in  Figure  55b,  in  which  the  bias 
begins  to  significantly  affect  closed  loop  tracking  synchronization  after  t  w  11  sec  («  17.5 
code  cycles).  For  these  reasons,  analysis  on  the  eMRDLL  was  conducted  for  the  two  time 
intervals  of  0  to  10  sec  and  0  to  30  sec;  simulations  were  not  conducted  for  longer  periods 
of  time  in  order  minimize  the  simulation  run  times. 


a)  a  =  0.3 


b)  a  =  0.5 


Figme  55  Discriminator  outputs  for  a  =  0.3  and  a  —  0.5. 


4.6.1  MCEU  Performance  During  Code  Phase  Tracking.  MCEU  estimates  of 
d,  ®0)  and  ®l  for  “  =  0-3  and  a  =  0.5  are  illustrated  in  Figures  56  through  57.  The 
figures  indicate  that  the  estimates  degreide  in  proportion  to  the  corresponding  discriminator 
output;  correspondingly,  the  degradation  is  less  significant  for  a  =  0.3  and  more  severe 
for  a  =  0.5.  The  eMRDLL  is  able  to  correct  itself  and  maintain  track  eis  indicated  by 
the  cyclic  nature  depicted  in  Figure  55.  This  cyclic  property  is  also  present  in  Figure  57. 
Although  a  bias  exists,  the  eMRDLL  does  not  loose  lock  as  a  result  of  this  bias. 


^Because  the  Simtilink  models  s^e  discrete  in  nature,  a  certain  VCC  input  threshold  must  be  exceeded 
prior  to  any  VCC  frequency  adjustment. 
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Figure  56  MCEU  estimates  a,  ®o>  ®l  for  a  =  0.3. 


Figure  57  MCEU  estimates  a,  xq,  and  ®i  for  a  =  0.5. 
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J^.6.2  eMRDLL  Versus  NCDLL  (A  =  1.0  and  A  =  0.1)  Noise-Free  Code  Phase 
Tracking.  Code  phase  error  tracking  measurements  were  taken  after  the  models  reached 
steady-state.  The  simulations  for  the  NCDLL  were  conducted  for  50  sec  to  ensure  that  at 
least  10,000  samples  were  available  for  determining  tracking  performance.  Shorter  simu¬ 
lation  nm  times  were  acceptable  for  the  eMRDLL,  because  simulations  were  started  with 
the  eMRDLL  initially  synchronized. 

Figure  58  provides  the  mean  code  phase  tracking  error  of  the  NCDLL  versus  the 
eMRDLL.  Because  the  operation  of  the  MCEU  is  constrained  to  estimating  multipath 
parameters  for  delays  a  G  [0.5, 1.5],  the  case  in  which  only  the  LOS  signal  exists  is  not 
considered  in  the  evaluation  of  the  eMRDLL.  The  envelopes  indicate  the  expected  theo¬ 
retical  worst-case  tracking  error  as  determined  by  the  equations  in  Table  1.  The  tracking 
performance  of  the  NCDLL  (A  =  1.0)  is  anedyzed  for  uniformly  spaced  delays  through 
the  range  of  of  G  [0,1.5]  in  increments  of  0.1;  tracking  performance  was  also  sinalyzed  for 
delays  a  =  0.55  and  a  =  0.75.  A  negative  tracking  error  of  the  NCDLL  corresponds  to  a 
reflected  signal  being  out  of  phase  with  the  LOS  signal,  as  is  the  case  here  for  a  —  0.55 
and  a  =  0.75.  As  expected,  the  NCDLL  with  correlator  spacing  A  =  0.1  outperforms 
the  NCDLL  with  correlator  spacing  A  =  1.0;  however,  the  results  of  the  NCDLL  with 
A  =  0.1  are  higher  {S  «  0.035)  through  a  G  [0.1, 1.0]  than  predicted  by  the  equations 
in  Table  1  (5  =  0.025).  Previous  anedysis  has  shown  that  a  NCDLL  with  a  normalized 
correlator  spacing  of  A  =  0.1  does  not  outperform  a  spacing  of  A  =  0.2  in  the  presence 
of  multipath  [18];  perhaps  the  failure  (of  the  NCDLL  with  A  =  0.1)  to  achieve  predicted 
performance  is  a  limitation  intrinsic  to  the  NCDLL.  As  illustrated  in  Figure  58,  the  mean 
code  phase  error  of  eMRDLL  is  lower  than  the  NCDLL  A  =  0.1  throughout  the  range  of 
a  G  [0.1, 1.0]  with  the  exception  of  a  =  0.5  (for  a  flnal  time  tf  —  10  sec).  For  tf  =  30  sec, 
the  code  phase  error  for  a  =0.5,  0.6,  and  0.7  is  worse  than  results  obtained  for  the  NCDLL 
(A  =  0.1),  indicating  that,  for  longer  periods  of  code  phase  tracking,  the  accumulating 
effect  of  the  bias  leads  to  greater  degradation  in  eMRDLL  tracking  performance. 

The  code  phase  tracking  rmse®,  commonly  known  as  tracking  jitter,  for  the  eMRDLL 
and  the  NCDLL  is  illustrated  in  Figure  59.  The  code  phase  tracking  rmse  for  the  NCDLL 

^The  rmse  reflects  the  efiects  of  the  bias  in  addition  to  the  variance  of  the  code  phase  tracking  error. 
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NCDLL  vs.  eMRDLL 


Figure  58  Code  phase  tracking  error  (mean)  where  represents  NCDLL  (A  =  1.0), 
‘+’  represents  NCDLL  (A  =  0.1),  ‘o’  represents  eMRDLL  (0  to  10  sec),  ‘x’ 
represents  eMRDLL  (0  to  30  sec). 

(A  =  1.0  and  A  =  0.1)  is  approximately  the  Scune  as  the  absolute  value  of  the  mean 
phase  error  which  indicates  that  negligible  bias  is  present  in  their  measurements.  The 
rmse  phase  error  of  the  eMRDLL  is  greater  than  the  mean  phase  error,  as  a  result  of 
the  bias.  Results  indicate  that  for  tf  =  10  sec,  the  eMRDLL  outperforms  the  NCDLL 
(A  =  0.1)  for  a  €  [0.1 ;  1.0]  with  the  exception  of  a  =  0.5.  As  with  the  mean  phase  error, 
the  eMRDLL  performance  degrades  when  tf  =  0  :30  sec  for  a  =  0.5,  0.6,  and  0.7  where  it 
is  outperformed  by  the  NCDLL  (A  =  0.1).  The  code  phase  tracking  rmse  for  other  values 
of  a  increase  only  slightly  due  to  the  tf  increase.  It  should  be  noted  that  the  NCDLL 
(A  =  0.1)  marginally  outperforms  eMRDLL  for  a  €  [1.1, 1.4);  this  is  a  reflection  of  the 
sampling  error  of  the  MCEU  as  demonstrated  in  Simulation  1  for  a  =  1.3  (in  which 
samples  were  taken  from  the  lower  boimd  of  the  MCEU  output).  This  sampling  error  is 
present  for  all  a  G  [1.1, 1.4].  The  eMRDLL  could  possibly  be  improved  by  reducing  its 
normalized  correlator  spacing  A.  Overall,  these  noise-free  baseline  results  indicate  that 
the  eMRDLL  performed  much  better  than  the  NCDLL  (A  =  1.0)  for  all  a  €  (0,1.5].  One 
may  speculate  that  if  Simulink  computational  error  is  solely  responsible  for  the  observed 
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bias,  then  the  eMRDLL  could,  uniformly  in  a,  outperform  the  NCDLL  (A  =  0.1)  in  an 
actual  receiver. 


NCDLL  vs.  eMRDLL 


Figure  59  Code  phase  tracking  rmse  where  represents  NCDLL,  A  =  1.0.  *+’  repre¬ 
sents  NCDLL,  A  =  0.1.  ‘o’  represents  eMRDLL  (0  to  10  sec),  ‘x’  represents 
eMRDLL  (0  to  30  sec). 


4.7  Simulation  #  4'  eMRDLL  Performance  in  AWGN 

Simulations  were  conducted  to  characterize  the  performance  of  the  eMRDLL  in  the 
presence  of  AWGN.  The  objective  was  to  determine  the  effects  of  different  SNR  levels  on 
the  eMRDLL  code  phase  tracking.  Five  reflected  signal  delays,  a  =  0.3,  0.5,  0.7,  1.0,  and 
1.3,  (the  same  delays  analyzed  in  Simulation  #  2)  were  used  for  simulation.  The  values  of 
a  chosen  for  analysis  represent  the  best  (a  =  0.3)  through  the  worst  (a  =  0.5)  code  phase 
tracking  results  of  Simulation  #  3.  As  before,  the  range  tf  :=  10  sec  was  used  for  analysis; 
however,  the  tf  :=  20  sec  was  used  instead  of  t/  :=  30  sec,  to  decrease  simulation  run 
times.  The  simulations  were  conducted  for  a  sampling  of  SNR€  [—35, 30]  dB,  as  indicated 
in  Table  5. 

The  mean  of  the  code  phase  tracking  error  for  both  ranges  of  t  is  provided  in  Fig¬ 
ure  60.  The  figure  illustrates  the  distinct  differences  in  the  fiuctuation  of  the  mean  errors 
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0.08 


10  second  interval 


20  second  interval 

1 - 1 - 1 - r 


Figure  60  Mean  of  code  phaise  tracking  error  for  a  =  0.3  (‘ — ’),  a  =  0.5  a  =  0.7 

(‘-  a  =  1.0  (‘.’)>  a  =  1.3 


over  the  different  intervals  of  time.  For  tf  :=  10  sec,  the  mean  error  is  relatively  constant 
for  all  delays  through  the  range  SNR6  [-13,30]  dB,  with  the  exception  of  a  =  0.5  (which 
was  shown  to  be  affected  by  the  bias  after  five  seconds;  see  Figure  55b).  At  SNR=-13  dB, 
the  mean  code  phase  tracking  error  for  a  =  1.0  begins  to  grow  as  a  result  of  the  noise. 
At  SNR=-35  dB,  the  mean  treicking  error  at  all  of  the  delays  has  been  affected  by  the 
noise.  For  tf  :=  20  sec,  the  mean  tracking  error  versus  SNR  fluctuates  more  dramatically 
for  a  =  0.5  and  a  =  0.7;  this  is  consistent  with  Figure  59,  wherein  a  =  0.7  was  shown 
to  have  a  significantly  larger  rmse  for  longer  simulation  runs.  The  mean  tracking  error 
of  the  remaining  delays  are  relatively  constant  throughout  the  range  SNRe  [—13,30]  dB. 
Note  that  for  a  =  0.3  the  tracking  error  did  not  show  degradation  until  SNR=-23  dB.  For 
a  =  0.5  and  a  =  0.7,  SNRG  [-13,-6]  dB  and  SNRG  [0,5]  dB,  respectively,  the  mean 
tracking  error  was  minimal;  this  illustrates  how  the  noise  exciting  the  system  can  actually 
mitigate  the  bias  effect.® 

Similar  results  were  observed  for  the  variance  of  the  code  pheise  tracking  error  as 
illustrated  in  Figure  61.  For  tf  :=  10  sec.  Figure  61a,  the  variance  cleeurly  exhibits  a 
SNR  threshold  of  -23  dB.  For  t  :=  20  sec,  this  is  also  true  with  the  exception  of  a:  =  0.7 
which  shows  an  erratic  variance  versus  SNR.  As  with  the  mean  error,  the  veiriance  of  the 

^This  is  known  as  a  ‘persistence  of  excitation’  condition. 
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a)  t=0:10  sec  b)  t=0:20  sec 

Figure  61  Variance  of  code  phase  error  for  a  =  0.3  (‘ — ’)>  “  =  0*5  (‘■*’))  ®  =  0*7  (‘-  -’)> 
a  =  1.0  and  a  =  1.3  (*-.’)• 


code  phase  tracking  error  for  a  =  0.5  and  a  =  0.7  was  actually  reduced  in  the  range  of 
SNRe  [0,30]  dB]  the  presence  of  noise  counteracts  the  detrimental  effects  of  the  bias  in 
this  range. 


Figure  62  Code  tracking  rmse  for  a  =  0.3  (‘ — ’),  a  =  0.5  ®  =  0-7  (‘-  -’)>  ®  =  1-0 

and  a  =  1.3 


An  alternative  measure  of  code  tracking  performance  is  the  code  phase  tracking  rmse 
or  tracking  jitter,  which  reflects  the  bias  and  Arariance  of  the  error  signal.  Figure  62  provides 
the  code  phase  rmse  for  tf  :=  10  sec  euid  tf  :=  20  sec.  For  tf  :=  10  sec,  the  rmse  of  the 
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code  phase  tracking  error  exhibits  a  threshold  at  SNR«  —23  dB.  For  tf  :=  20  sec,  the  SNR 
performance  threshold  drops  to  -13  dB  with  the  exception  of  a  =  0.5  and  a  =  0.7  which 
display  erratic  rmse  behavior  versus  SNR. 

In  summary,  the  ctddition  of  moderate  levels  of  noise  can  cictually  improve  the  per¬ 
formance  of  the  eMRDLL.  The  noise  has  the  effect  of  perturbing  or  exciting  the  system 
such  that  the  effects  of  the  Simulink  bias  are  reduced.  Of  course,  there  is  a  noise  threshold 
for  which  the  code  phase  tracking  of  eMRDLL  begins  to  degrade.  Prom  the  figures  in  this 
section,  one  can  see  that  SNR  level  at  which  the  performance  of  the  eMRDLL  degrades 
varies  for  different  delays;  however,  generally  an  SNR  of  -23  dB  (relative  to  the  Simulink 
models)  is  the  approximate  SNR  level  at  which  significantly  degradaded  code  tracking  per¬ 
formance  can  be  expected  for  tf  :=  10  sec.  An  SNR  threshold  of  -13  dB  (relative  to  the 
Simulink  models)  is  expected  when  the  interval  for  averaging  is  doubled. 

As  explained  in  the  summary  of  Simulation  4^  2,  the  specified  SNR  level  that  would 
be  equivalent  to  GPS  SNR  levels  is  »  —35  dB.  Results  show  that  at  this  SNR,  significant 
degredation  in  tracking  performance  occurs.  Reducing  the  bandwidth  of  the  LPFs  and 
increasing  the  length  of  the  PN  sequence  to  N=1023  may  improve  performance. 
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J^.8  Simulation  #  5:  Time  Varying  Reflections 

Previous  simulations  were  conducted  under  the  assumption  that  the  delay  of  the 
reflected  signal  was  slowly  time  varying,  as  would  be  encountered  with  static  position  de¬ 
terminations.  Under  this  condition,  it  was  assumed  that  the  multipath  delay  over  relatively 
short  periods  of  time  was  constant. 

Many  realistic  conditions  involve  time  varying  multipath  delays.  Simulations  were 
conducted  on  the  MCEU  given  a  time  varying  multipath  environment,  wherein,  it  was 
assumed  that  the  MCTL  maintained  perfect  synchronization  with  the  received  signal.  In 
addition,  it  was  assumed  that  the  received  signal  has  already  been  perfectly  filtered  to 
baseband  (i.e.,  carrier  modulation  was  removed). 

Two  time  varying  multipath  scenarios  aire  considered. 

•  Case  1:  a{t)  —  (10/63)t  —  2.17/63  for  t  =  0  : 10  sec 

•  Case  2:  a(t)  =  (5/63)t  -  1.085/63  for  t  =  0  :  19.5  sec 

Time  duration  is  chosen  to  ensure  that  delays  of  a  €  [0, 1.5]  are  encountered  within 
the  simulation  time.  Estimates  are  calculated  by  the  MCEU  every  0.63  sec. 

Figure  63  illustrate  the  results  of  the  simulations  conducted  for  case  1.  Figure  63a 
indicates  that  the  estimate  of  a  is  lower  than  the  true  for  all  a,  with  the  exception  of 
a  =  0.1.  This  under  estimation  of  a  is  likely  a  result  of  the  narrow  bandwidth  of  the  LPFs 
in  the  bank  of  correlators,  because  the  narrow  bandwidth  LPFs  invoke  a  time  delay  on  the 
filtered  signal;  consequently,  current  estimates  are  based  upon  previous  correlator  samples 
which  have  been  delayed  by  the  lag/response  of  the  LPF.  The  figure  also  indicates  that 
XQ  and  x\  are  very  accurate  with  the  exception  of  the  estimates  taken  during  the  first  two 
cycles,  which  is  the  transient  period.  Note  that  the  increased  error  in  the  estimates  xq  and 
®1  for  Q!  €  [1.1, 1.5]  is  consistent  with  results  presented  in  Simulation  #  1. 

Because  of  the  lag  in  the  MCEU  estimates  observed  in  the  first  case,  a  slower  time 
varying  simulation  was  chosen  for  case  2,  which  has  a  delay  which  increases  at  half  the 
rate  of  the  first  case.  As  illustrated  in  Figure  64,  with  the  exception  of  the  first  estimate, 
the  estimates  of  a  are  indeed  closer  to  the  true  delays,  as  expected,  since  the  change  in  a 
over  time  is  less  rapid.  Note  that  the  estimates  xq  and  xi  are  nearly  identical  to  those  of 
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Figure  63  Time  varying  multipath  estimates  a,  xq,  and  xi. 


the  first  case,  which  suggests  that  the  estimates  in  the  range  of  a  6  [1.1, 1.5]  are  degraded 
due  to  the  instantaneous  value  of  alpha  as  opposed  to  any  time  varying  properties  of  the 
reflected  signal. 


Figure  64  Time  varying  multipath  estimates  a,  ®o»  and  xi. 


The  effects  of  noise  upon  time  varying  multipath  estimation  were  also  considered. 
The  SNR  levels  considered,  0  dB  and  —23  dB,  were  selected  because  previous  results 
indicate  that  slight  noise  degradation  occurs  at  0  dB  and  significant  degradation  occurs 
at  —23  dB  (the  SNR  threshold).  Figure  65  illustrates  that,  as  expected,  a  0  dB  SNR  has 
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little  effect  on  the  estimates  of  a.  However,  there  is  a  slight  degradation  in  the  estimates 
XQ  and  x\. 


Figure  65  Time  varying  multipath  estimates  d,  ®o>  aJid  ®i  for  SNR=0  dB. 


Figure  66  illustrates  that  decreasing  the  SNR  level  to  -23  dB  does  begin  to  signifi¬ 
cantly  effect  the  estimates  of  a,  xq  and  xi,  as  expected.  The  overall  results  indicate  that 
the  performance  of  the  MCEU  with  time  varying  multipath  is  consistent  with  the  perfor¬ 
mance  encountered  under  a  constant  multipath  delay,  both  with  and  without  the  presence 
of  AWGN. 


Figure  66  Time  varying  multipath  estimates  a,  xq,  etnd  ®i  for  SNR=0  dB. 
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V.  Conclusion  and  Recommendations 


5.1  Overview 

This  thesis  proposed  the  use  of  maximum  likelihood  estimation  in  a  code  traxdcing 
loop  design  to  mitigate  the  effects  of  multipath.  Although  the  NCDLL  with  a  normalized 
narrow  correlator  spacing  of  A  =  0.1  has  been  showed  to  reduce  code  phase  tracking  error 
caused  by  multipath  relative  to  a  standard  (A  =  0.1)  NCDLL,  it  cannot  fully  remove  the 
tracking  error.  The  eMRDLL,  which  is  an  enhanced  version  of  the  MRDLL  [6,7],  accounts 
for  the  parameters  of  the  reflected  as  well  as  the  direct  path  signals  to  mitigate  the  effects 
of  multipath.  The  effectiveness  of  the  eMRDLL  is  highly  dependent  on  the  performance 
of  the  MCEU,  which  is  a  MVU  estimator  capable  of  reaching  the  CRLB. 

Analysis  conducted  on  the  equations  which  describe  the  operation  of  a  standard 
NCDLL  characterized  a  tracking  curve  known  as  the  S-curve.  Given  a  direct  path  only 
signal,  the  steady  state  tracking  point  of  this  S-curve  is  located  at  8*  =  0;  however,  in 
a  multipath  environment,  the  S-curve  becomes  distorted  and  treicks  with  a  bias  6*  ^  0. 
Further  analysis  confirmed  that  reducing  the  correlator  spacing  from  A  =  1.0  to  A  =  0.1 
reduced  the  effects  of  multipath  significantly;  however,  a  threshold  was  reached  for  which 
tracking  errors  could  not  be  reduced  further  due  to  inherent  constraints  of  the  NCDLL 
design. 

The  MRDLL,  which  is  composed  of  the  MCTL,  MCEU,  and  the  ALC,  has  previously 
been  shown  to  mitigate  the  effects  of  multipath  [6,7];  however,  operation  of  the  MRDLL  was 
limited  to  estimating  multipath  delays  through  the  range  of  a  €  [0.1, 1,5]  in  increments  of 
0.1.  The  eMRDLL  provides  enhancements  to  the  MRDLL,  via  an  enhcinced  MCEU,  which 
permits  multipath  parameters  to  be  estimated  through  the  full  range  of  a  6  [0.1, 1.5].  An 
analysis  was  conducted  to  examine  the  performance  improvements  of  the  eMRDLL  in  a 
multipath  environment.  Anedysis  of  the  MCTL  loop  equations  indicate  that  a  tracking 
curve  (S-curve)  is  formed  from  a  combination  of  operating  curves  based  on  the  direct  path 
and  reflected  path  signal  parameters,  the  sum  of  which  generates  a  steady-state  tracking 
point  of  =  0  when  the  MCEU  estimates  axe  exact. 
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Analysis  was  conducted  on  the  enhanced  MCEU  in  order  to  characterize  its  perfor¬ 
mance.  A  SQP  search  algorithm  locates  the  minimum  of  the  quEidratic  cost  function  V{a), 
maximizing  the  PDF  of  the  linear  statistical  measurement  model. 

5.2  Computer  Simulations 

Simulations  were  conducted  using  Simulink  to  determine  the  tracking  performance 
of  the  eMRDLL  in  a  realistic  multipath  environment.  Five  separate  simulations  were 
conducted.  The  first  simulation  focused  on  the  stand  alone  noise-free  performance  of  the 
MCEU  (assumes  perfect  MCTL  synchronization).  The  second  simulation  was  conducted  to 
analyze  the  noise-free  closed  loop  code  phase  tracking  performance  of  the  eMRDLL  versus 
the  tracking  performance  of  the  NCDLL  with  norm«dized  correlator  spacings  of  A  =  1.0 
and  A  =  0.1.  The  third  simulation  focused  on  the  stand  alone  performance  of  the  MCEU 
in  the  presence  of  AWGN  (assumes  perfect  MCTL  synchronization).  The  fourth  simulation 
was  conducted  to  analyze  the  closed  loop  code  phase  trricking  performance  of  the  eMRDLL 
in  the  presence  of  AWGN.  Finally,  simulations  were  conducted  on  the  stand  alone  MCEU 
(assumes  perfect  MCTL  synchronization)  in  a  time  varying  multipath  environment,  with 
and  without  AWGN. 

Note:  SNR  performance  is  relative  to  the  Simulink  environment  (see  Section  4.3.7). 

5.2.1  Simulation  #  1  Results.  Simulations  of  the  noise-free  stand  alone  MCEU 
resulted  in  accurate  estimates  of  the  multipath  parameters  within  a  multipath  environ¬ 
ment.  The  MCEU’s  performance  is  sensitive  to  the  LPFs  used  in  the  bank  of  correlators. 
The  LPFs  function  as  integrators,  convert  the  received  signal  to  baseband  (an  enhancement 
provided  by  the  eMRDLL),  and  attenuate  noise.  The  response  of  the  filters  implemented 
contributed  significantly  to  the  eiccureicy  of  the  estimates  provided  by  the  ML  SQP  algo¬ 
rithm. 

5.2.2  Simulation  #  2  Results.  Simulations  of  the  stemd  alone  MCEU  in  the 
presence  of  AWGN  indicate  that  it  performs  well.  For  the  models  tested,  decreased  per¬ 
formance  in  estimation  was  not  evident  until  the  SNR  was  reduced  to  0  dB.  Significant 
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degradation  in  MCEU  performance  was  not  realized  until  the  SNR  was  reduced  to  the 
identified  threshold  of  —23  dB. 

5.2.3  Simulation  #  3  Results.  Simulations  indicate  that  a  bieis  is  generated  which 
progressively  hinders  the  performance  of  the  MCTL  over  time.  Although  the  bias  is  present 
(likely  due  to  nonidealities  of  the  Simulink  implementation),  code  phase  tracking  of  the 
eMRDLL  is  still  typically  better  than  code  phase  tracking  of  the  baseline  narrow  correlator 
NCDLL  with  for  a  G  [0.1, 1.0],  with  the  exception  of  delays  a  =  0.5,  a  =  0.6,  and  a  =  0.7. 
It  is  speculated  that  successful  removal  of  the  bias  would  enhance  performance  of  the 
eMRDLL  such  that  it  would  outperform  the  NCDLL  (A  =  0.1)  for  all  delays  a  G  [0.1, 1.0]. 
The  narrow  correlator  NCDLL  exhibited  slightly  better  performance  than  the  eMRDLL 
for  delays  a  6  [1.1, 1.5];  consequently,  code  tracking  performemce  of  the  eMRDLL  for  these 
delays  could  possibly  be  improved  by  reducing  the  correlator  spacing  within  the  MCTL  to 
A  =  0.1.  The  eMRDLL  outperforms  the  NCDLL  with  correlator  spacing  A  =  1.0  for  all 
delays  a  G  [0.1, 1.5]. 

5.2..^  Simulation  #  4  Results.  Simulations  of  the  eMRDLL  in  the  presence  of 
AWGN  indicate  that,  as  with  the  MCEU,  significant  degradation  in  code  tracking  perfor¬ 
mance  does  not  occur  until  the  SNR  is  reduced  to  —23  dB.  Furthermore,  small  levels  of 
noise,  SNR  «  0  dB,  enhanced  the  performance  of  the  eMRDLL  relative  to  the  noiseless 
case  by  exciting/perturbing  the  system  to  counter  the  effects  of  the  bias.  The  tracking 
quality  had  a  tendency  to  degrade  over  time,  due  to  the  cummulative  effect  of  the  bias. 

5.2.5  Simulation  #  5  Results.  Simulations  conducted  on  the  MCEU  in  a  time 
varying  multipath  environment  indicate  that  the  MCEU  provides  estimates  which  are  very 
accurate  and  comparable  with  those  for  the  cases  of  fixed  multipath  delays.  There  is 
typically  some  estimation  error  introduced  by  the  signal  delay  caused  by  the  response  of 
the  LPF.  As  with  previous  simulations,  the  MCEU  showed  no  degredation  in  estimates 
for  SNR>  0  dB.  Significant  degredation  in  multipath  parameter  estimation  occurred  for 
SNR<  -23  dB. 
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5.3  Summary  and  Recommendations 

This  thesis  has  shown  that  the  eMRDLL  has  tremendous  potential  for  multipath 
mitigation  in  GPS  applications.  Although  the  parameters  used  in  the  simulations  were 
different  than  those  found  in  GPS,  the  code  tracking  loop  models  exhibited  properties  that 
were  very  similar  to  the  properties  encountered  in  the  GPS  environment.  The  eMRDLL 
exhibited  significant  improvements  in  steady  state  code  phase  tracking  in  comparison  with 
the  code  tracking  performance  of  the  standard  NCDLL  in  the  presence  of  a  single  refiected 
signal.  Furthermore,  the  eMRDLL  exhibited  improved  code  phase  tracking  in  comparison 
with  the  code  tracking  performance  of  the  narrow  correlator  NCDLL  for  most  delays 
a  €  [0.1, 1.0].  Removal  of  the  simulation  bias  generated  by  the  Simulink  implementation 
should  further  increase  performance.  The  MLE  implemented  in  the  eMRDLL  has  been 
shown  to  accurately  estimate  the  multipath  signal  parameters  for  delays  a  €  [0.2, 1.5]. 
The  following  is  a  list  of  recommendations  for  follow-on  research  to  further  improve  the 
performance  of  the  eMRDLL. 

•  Reduce  the  correlator  spacing  within  the  eMRDLL  to  improve  code  phase  tracking 
performance  for  delays  a  €  [1.05, 1.5]. 

•  Generalize  the  eMRDLL  for  operation  in  a  mvdtipath  environment  where  multiple 
reflections  are  present. 

•  Examine  the  effects  of  Doppler  on  the  code  phase  of  the  received  multipath  signal 
with  and  without  the  ALC  implemented. 

•  Examine  the  performance  of  the  eMRDLL  when  random  data  modulation  is  in¬ 
cluded  on  the  received  GPS  signal,  and  perform  bit  error  analysis. 

•  Design  a  detection/estimation  scheme  that  detects  the  presence  of  multipath.  This 
would  allow  the  NCDLL  to  be  used  when  multipath  is  not  present,  since  it  has  been  shown 
to  perform  better  in  the  absence  of  multipath  [6] . 
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Appendix  A.  Computer  Simulation  Models 

A.l  Overview 

Two  basic  models  were  designed  using  Simulink:  the  eMRDLL  and  the  NCDLL. 
Parameters  used  for  the  simulations  are  provided  in  Section  4.3.8.  This  appendix  provides 
a  description  of  the  Simulink  models  designed. 

A. 2  The  Transmitted  Multipath  Signal 

The  multipath  environment  was  simulated  by  creating  a  PN  code  generator  which 
operated  at  a  rate  governed  by  a  pulse  generator  (l/rc=100  Hz)  as  shown  in  Figure  67. 


Figure  67  Simulink  multipath  model. 

The  PN  generator  consists  of  a  six  bit  feedback  register,  shown  in  Figure  68,  with  a 
polynomial  X^+X^+X^+X^+1  of  order  m  =  Q.  This  register  produces  a  maximum  lenth 
PN  sequence  of  iV  =  2"*  —  1  =  63  [3].  Although  GPS  transmits  a  signal  with  code  length 
N  =  1023,  choosing  N  =  63  for  the  simulations  still  permitted  the  use  of  the  large  N  code 
correlation  approximation  of  Equation  8  with  the  added  benefit  of  increasing  simulation 
performance. 

Eanh  register  of  the  PN  generator,  shown  in  Figme  69,  consists  of  switches  emd  delays, 
which  are  controlled  by  the  pulse  rate  of  the  PN  generator.  Each  register  is  initiated  with 
initial  setting  equal  to  one.  This  is  accomplished  by  setting  the  initial  condition  of  the  unit 
delay  block  in  eanh  register  equal  to  one.  The  received  pulses  from  the  pulse  generator 
arrive  at  input  2,  which  flips  the  switch  allowing  the  next  bit  to  enter  through  input  1. 
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Figure  68  Simulink  PN  generator  model. 

This  value  then  flips  switch  1  to  the  corresponding  constant  value  which  departs  through 
output  1  and  also  cycles  through  the  loop  until  the  next  pulse  flips  the  switch  again.  As 
such,  the  output  of  each  register  cycles  through  the  next  register  upon  receiving  a  pulse 
from  the  pulse  generator. 


Figure  69  Simulink  model  of  a  single  register,  R,  of  the  PN  code  generator. 

The  output  chip  from  the  PN  generator  then  proceeds  to  the  switch  shown  in  Fig¬ 
ure  68  which  flips  the  switch  to  -1-1  for  code  value  equal  to  1  or  -1  for  code  value  equal 
to  0.  Thus  generating  a  sequence  of  ±ls.  This  sequence  is  then  modulated  onto  a  carrier 
/o=1000  Hz  with  phase  equed  to  7r/2  (cosine  wave).  The  sequence  of  ±ls  are  represented  as 
180°  phase  shifts  in  the  carrier  signal.  This  represents  the  transmitted  direct  path  signal. 
The  reflected  signal  is  then  generated  by  delaying  and  attenuating  the  direct  path  signal 
and  then  summing  this  signal  with  the  direct  path  signal  to  form  the  multipath  signal. 
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A.2.1  NCDLL  Model.  The  simulink  model  of  the  NCDLL  is  beised  on  the  NCDLL 
illustrated  in  Figure  4.  Figure  70  shows  the  model  implemented  in  Simulink.  The  figure 
indicates  that  the  transmitted  signal  is  delayed  prior  to  entering  the  ‘early-late’  gate  chan¬ 
nels  of  the  NCDLL.  This  delays  the  transmitted  signal  relative  to  the  ‘early’  signal  of  the 
NCDLL  by  ATc/2  and  mainteiins  an  advance  of  the  transmitted  signal  relative  to  the  ‘late’ 
signal  by  ATcf2,  The  remaining  operations  of  the  NCDLL  are  based  on  the  theory  of 
operations  described  by  the  equations  in  Chapter  2  for  the  NCDLL.  Parameters  chosen  for 
the  model  are  listed  in  Section  4.3.8. 


Figure  70  Simulink  model  of  the  NCDLL. 


A.S.2  Simulink  Models  of  the  VCC  and  the  Local  PN  Generator.  The  VCC 
implemented  in  the  NCDLL  is  shown  in  Figure  71.  The  error  signal  from  the  output 
of  the  loop  filter  enters  the  VCC  and  is  multiplied  by  a  gain  of  -1.  This  is  necessary 
in  order  for  proper  synchronization  to  occur.  The  signal  then  enters  the  VCO  (block 
provided  by  the  Communications  Toolbox)  which  operates  a  given  quiescent  frequency 
(100  Sz).  The  frequency  of  the  VCO  is  adjusted  according  to  the  input  error  signal. 
The  sinusoidal  output  of  the  VCO  then  enters  a  rising  edge  detector  (block  provided  by 
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the  Communications  Toolbox)  which  produces  a  pulse  at  the  first  positive  leading  edge 
of  the  output  of  the  VCO.  The  rate  of  these  pulses  determines  the  code  rate  of  the  local 
PN  generator  shown  in  Figure  72.  Note  that  the  local  PN  generator  operates  the  same 
way  as  the  transmitted  signal  generator  previously  described.  If  the  received  signed  is 
delayed  (out  of  synchronization)  with  the  locally  generated  signal  then  the  discriminator 
output  will  be  some  positive  voltage  level  as  described  by  the  S-curve  for  the  NCDLL. 
After  passing  through  the  loop  filter,  this  signal  is  multiplied  by  a  negative  gain  of  one  and 
passed  through  to  the  VCO.  The  negative  signal  input  to  the  VCO  reduces  the  frequency 
of  the  sinusoidal  outputs  of  the  VCO  which  reduces  the  rate  of  the  clock  pulses  produced 
by  the  rising  edge  detector.  As  such,  the  rate  of  the  local  PN  generator  is  reduced  to  match 
the  delayed  received  signal. 


Figure  71  Simulink  model  of  the  VCC. 


Figure  72  Simulink  local  PN  generator  model. 
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A.S  eMRDLL 


The  simulink  model  of  the  eMRDLL,  illustrated  in  Figure  73  is  based  on  the  eMRDLL 
illustrated  in  Figure  15.  The  basic  components  of  the  Simulink  model  include  the 

•  transmitted  signal  and  transmission  delay  (previously  discussed), 

•  AWGN  function  block  (block  provided  by  the  Communications  Toolbox), 

•  PLL  model, 

•  direct  and  reflect  path  ‘eaxly-late’  gate  models, 

•  ALC  model, 

•  MCEU  model, 

•  •variable  delay  function  block, 

•  VCC  (pre'viously  discussed), 

•  local  PN  generator  model  (previoulsy  discussed). 


Figure  73  Simulink  model  of  the  eMRDLL. 


A.3.1  Simulink  AWGN  Block.  The  AWGN  block  is  a  Simulink  function  block 
provided  by  the  Communications  Toolbox.  There  are  three  parameter  settings  for  this 
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block:  the  mean,  the  seed  value,  and  the  variance.  The  mean  weis  always  set  to  zero.  The 
seeds  were  varied  for  different  SNRs  as  specified  in  Table  7.  The  variance  of  the  noise  was 
used  in  determining  the  SNR  as  expledned  in  Section  4.3.7. 

Table  7  Specified  SNR  with  corresponding  seed  value. 


Specified  SNR  (dB) 

Seed  Value  (dB) 

30 

12345 

25 

12345 

20 

12345 

15 

46531796 

10 

12345 

5 

46531796 

0 

46531796 

-6 

84531796 

-13 

84531796 

-23 

84531796 

-35 

84531796 

A.S.2  Simulink  PLL  Model.  The  received  signal  is  mixed  with  a  signal  from  the 
PLL  which  is  designed  as  a  black  box  that  provides  the  theoretical  PLL  signal  as  described 
in  Chapter  3.  The  mixed  signal  then  enters  the  MCEU  and  the  ‘early-late’  gates  of  the 
MCTL. 

A.3.3  Simulink  ‘Early-Late’  Gate  Model.  The  ‘early-late’  gate  design  of  the 
direct  and  refiected  paths  is  illustrated  in  Figure  74.  This  is  a  typical  coherent  design  in 
which  the  received  signal  is  converted  to  baseband  and  mixed  with  a  locally  generated 
code  replica  spaced  at  ±ATc/2  in  the  ‘early’  and  ‘late’  channels.  Low  pass  filtering  the 
signals  and  differencing  provides  the  discriminator  in  each  channel  as  defined  in  Chapter  3. 
The  signals  pass  through  the  gain-phase  correlators  and  are  summed  to  form  the  eMRDLL 
discriminator  output  (S-curve). 

A.3.4  Simulink  ALC  Model.  The  discriminator  output  signal  enters  the  ALC 
illustrated  in  Figure  75  which  is  based  on  the  Figure  22  and  Equations  92  and  93.  Ka  = 
1/ti  =  8/225A  is  determined  by  Equation  94  where  a;„  «  O.OGtt  rad/ sec,  Ko  =  1,  and 
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Figure  74  Simulink  ‘early-late’  gate  model  used  in  the  eMRDLL. 

If  A  is  determined  by  the  matlab  function  alc.m  based  on  the  inputs  from  the  MCEU. 
T2  —  2C/un  =  7.5  where  ^  =  l/V^  remziined  fixed  as  with  the  NCDLL.  The  output  of  the 
ALC  drives  the  VCC  which  drives  the  loccd  PN  generator  as  previously  described  for  the 
NCDLL. 


Figure  75  Simulink  model  of  the  ALC. 


A,3. 5  Simulink  MCEU  Model.  The  mixed  signal  entering  the  MCEU  is  converted 
to  baseband  and  correlated  with  locally  generated  PN  code  incrementally  spaced  by  BkTc. 
The  MCEU  implemented  in  Simulink  is  illustrated  in  Figure  76.  Two  inputs  exist  into 
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the  MCEU.  The  received  (mixed)  signal  enters  the  MCEU  through  input  1.  The  locally 
generated  spreading  code  enters  the  MCEU  through  input  2.  Note  that  the  locally  gen¬ 
erated  code  is  initially  delayed  by  the  same  delay  as  the  transmitted  delay.  This  ensures 
that  the  locally  generated  code  is  synchronized  with  the  transmitted  code.  The  locally 
generated  code  then  enters  a  bank  of  correlator  delays  which  are  incrementally  spaced  by 
0.15  through  the  range  Pk  6  [0,1.5].  the  delayed  code  then  enters  the  bank  of  correlators 
were  it  is  mixed/multiplied  with  the  received  (mixed)  signal  and  low  pass  filtered.  The 
LPFs  serve  to  convert  the  signsil  to  beiseband  and  also  operate  as  integrators.  The  corre- 


delay 

Figure  76  Simulink  model  of  the  MCEU. 

lator  samples,  R,  enter  the  MLE  sample  controller  block  shown  in  Figure  77.  The  sample 
controller  block  is  controlled  by  the  pulse  generator  which  is  connected  to  the  switch.  The 
pulse  generator  sends  a  pulse  at  the  sampling  rate  chosen  for  the  MLE,  NTc  =  0.63.  The 
switch,  when  triggered  by  the  pulse  generator,  allows  the  correlator  samples  R  to  pass 
through  at  the  appropriate  sampling  epoch.  Otherwise,  the  switch  passes  inputs  from  the 
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mux  which  contains  three  memory  blocks  for  a,  ®o>  These  memory  blocks  contain 

feedback  from  the  MLE.  The  general  concept  of  the  memory  block  is  to  pass  new  correlator 
samples  to  the  estimator  every  estimator  sampling  period,  NTc.  Otherwise,  the  previous 
estimates  of  a,  xq,  and  are  passed  through  the  MLE.  The  MLE,  which  is  written  as 
a  Matlab  function  is  designed  to  calculate  new  estimates  d,  io>  x\  when 

new  correlator  samples  are  received  and  pass  through  previous  unchanged  estimates,  which 
helps  to  reduce  processing  time.  The  MLE  calculates  estimates  by  using  an  SQP  algorithm 
as  described  in  Chapter  3. 


Figure  77  Simuhnk  model  of  the  sample  controller  block. 


Estimates  from  the  MCEU  axe  provided  to  three  different  functional  blocks  within 
the  eMRDLL:  the  ALC,  the  gain-phase  multipliers,  and  the  variable  delay  function  block. 


Figure  78  Simulink  model  of  the  variable  delay  function  block. 


^This  function  requires  the  Optimization  Toolbox. 
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A. 3.6  Simulink  Variable  Delay  Model.  The  variable  delay  function  block,  shown 
in  Figure  78,  controls  the  delay  of  spreading  code  from  the  local  PN  generator  to  ensure 
proper  synchronization  within  the  MCTL  reflected  path  channel.  The  locedly  generated 
spreading  code  enters  the  variable  delay  function  block  through  input  1  and  is  input  into 
two  multipath  variable  delay  blocks.  Each  of  the  two  multipath  veiriable  delay  blocks, 
shown  in  Figure  79,  is  comprised  of  the  same  components.  These  components  consist 
of  several  switches  of  which  any  combination  of  delay  a  G  (0,1.5]  in  increments  of  .01 
can  be  formed.  The  individual  delay  blocks  are  controlled  by  switches.  These  switches 
are  controlled  by  the  variable  thresholds  block  located  within  the  variable  delay  function 
block.  The  variable  thresholds  block  is  a  Matlab  function  (vdelay.m)  which  controls  the 
combination  of  switches  in  the  multipath  delay  block  beised  on  the  estimated  delay  a 
provided  by  the  MCEU. 


Figure  79  Simulink  model  of  the  multipath  variable  delay  block. 


The  output  of  the  two  multipath  delay  blocks  are  controlled  by  a  pulse  generator 
connected  to  a  switch  which  is  set  such  that  outputs  from  the  two  blocks  are  alternately 
passed  through  the  switch  every  other  estimator  sample  cycle  NTc.  Two  separate  multipath 
delay  blocks  are  used  to  ensure  that  all  delayed  chip  samples  within  the  combination  of 
delay  blocks  of  the  multipath  variable  delay  block  are  properly  cleared  when  a  new  estimate 
of  a  is  provided  by  the  MCEU. 
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A.J^  Matlah  Functions  Used  in  Simulink  Models 

The  following  Matlab  m-files  are  functions  used  in  the  Simulink  components  of  the 
eMRDLL. 

A.4’i  mpllm. 

function  signal^  mpll (parameters) 

*/,  MPLL  is  a  function  which  produces  a  signal  that  models  the  PLL  for  MRDLL. 

'/,  The  input  parameters  provides  t,  alpha,  aO,  euid  al. 

'/,  t  is  the  time  insteuit.  alpha  is  the  multipath  delay.  aO  is  the  amplitude  of 
'/,  the  direct  path  signal,  al  is  the  amplitude  of  the  reflected  signal. 

*/,  Authored  by  Fred  Baier 
'/,  Created  20  Oct  97 

t=parameters(l) ; 
alpha=parameters(2) ; 
aOsparameters (3) ; 
al=parameters (4) ; 

f_c=1000; 

T_c=l/100; 

theta0=0; 

if  alpha  <=  1 
R_c=l-alpha; 
else 
R_c=0; 
end 

A.4’2  alc.m. 
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function  Gain<salc  (estimates) 

*/,  ALC  is  a  function  that  determines  the  gain  of  the  Adaptive  Loop  Controller 
'/,  in  the  Simulink  environment  based  on  the  estimates  of  alpha,  xO,  and  xl. 

*/,  Authored  by  Fred  Baier 
'/,  Created:  17  Sep  97 
•/,  Modified:  9  Oct  97 

alpha=estimates(l) ; 
x0=estimates(2) ; 
xl=estimates(3) ; 

if  alpha  <*0.5 

A  =  2>i'(x0“2+xl~2+2*x0*xl); 
elseif  alpha  ==0.5 

A  =  2*(x0‘'2+xl~2+.5*x0*xl) ; 
elseif  alpha  ==1.6 

A  =  2*(x0“2+xl“2-.5*x0*xl) ; 
else 

A  =  2*(x0~2+xl''2-x0*xl) ; 
end 

Gain=l/A; 

A.4-3  mle.m. 

function  output=mle(R_in) 

'/function  [alpha , x_0 ,  x_l]  =mle  (R) 

'/.  MLE 

'/,  Authored  by  Fred  Baier 
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%  Functions  required:  crosscor.m,  toeplitz,  constr.m,  sigl.H.m,  V.alpha.m 


%  Function  created  for  simulink:  25  Aug  97 
%  Last  modified:  3  Sep  97,  22  Oct  97 

if  sum(R_in(4:ll))==80 

output>R_in(l:3) ;  '/,  NOTE:  Output  will  be  a  vector 

else 

*/.  THE  FOLLOWING  CALCULATES  NEW  OUTPUT  VECTOR 
global  R  Beta  C  M; 

R=R_in; 

7,  Create  noise  coveoriance  matrix,  C. 

M-11;  \7i  number  of  arms  in  MRDLL 

increment=0 . 15 ;  7.  delay  spacing  between  arms 

Beta= [0 : increment : 1 . 5] ’ ; 

Beta_diff=Beta; 

R_c“crosscor(Beta_dif f ,M) ; 

C=toeplitz(R_c) ; 

7i  STEP  1  is  performed  by  function  V.alpha. 

7>  global  R  Beta  C  M; 

7,  STEP  2.  Invoke  an  optimization  routine. 
f_old»le6:  %  Start  with  large  number  as  placeholder 
alpha_old=2 ; 
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for  2arin=l :  10 


if  €urm““l 

start_alpha=.l; 

VLLB=.05; 

VLUB=. 14999; 

•/.VLLB*.06; 

•/.VLUB=.149; 
elseif  aria“2 
steu:t_alpha= .  2 ; 

VLLB=. 15001; 

VLUB=. 29999; 

•/.VLLB-.161; 

•/,VLUB=.299; 

end 

junk=f opt ions; 

alpha=constr(’V_alpha’ ,start_alpha, junk,VLLB,VLUB) ; 

H_min=sigl_H (Beta, alpha, M) ; 

f = (inv (C) *R_in) ’ * (C-H_min*inv(H_min ’ *inv (C) *H_min) *H_min* ) * (inv(C) *R_in) ; 

if  f_old  <  f 

alpha*alpha_old ; 
else 

f _old=f ; 

alpha_old=alpha ; 
end 

’/.alpha  '/.  ADDED  TO  TRACK  IN  SIMULINK 

start_alpha=staxt_alpha+ . 16 ; 

VLLB=VLLB+.15; 
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VLUB=VLUB+.15; 

end 

R_in;  V.  FOR  TRACKING  IN  SIMULINK 
alpha; 

H=sigl_H(Beta, alpha, M) ; 

x_est*inv (H ’ *inv (C) *H) *H ’ *inv(C) *R_in ; 

x_0=x_est(l) ; 

x_l=x_est(2) ; 

output= [alpha  x_0  x_l] ’ ; 

end 


A.4-3.1  crosscor.m. 
function  R_c=crosscor (Omega, M) 

'/,  CROSSCOR  produces  cross  correlation  for  input  Omega 
for  index_a=l:M 

if  abs(Omega(index_a,l))  <=  1 

R_c ( index.a , 1 ) =l-abs (Omega ( index.a , 1) ) ; 
else 

R_ c ( index_a , 1 ) =0 ; 
end 
end 


A.4-3.2  siglJS.m. 
function  Hssigl_H(Beta, alpha, M) 

'/,  SIGl.H  produces  the  regressor  matrix,  H,  of  signal  for  vector  input  Beta,  (a 
'/,  vector  of  delays  for  each  arm),  alpha  (delay  of  the  reflected  signal,  single 
y.  reflection  case) ,  and  M  (the  number  of  arms  in  the  DLL) .  Requires 
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'/•  the  use  of  the  function  file  CROSSCOR. 
y,  Authored  by  Fred  Baier. 

Omega_l=Beta; 

R_cl=crosscor(Omega_l,M) ; 

0mega_2=alpha-Beta; 

R_c2=crosscor(0mega_2,M) ; 

H=[R_cl  R_c2]; 

A.4-3.3  V.alpha.m. 
function  Cf,g]=V_alpha( alpha) 

X  V.ALPHA  is  a  function  which  determines  alpha  by  minimizing  the  MLE  equ. 

'/•  Want  to  minimize  MLE  with  constraints  on  alpha. 

'/,  R  is  correlator  output  data  (known)  .  Beta  is  the  delay  of  the  arms  of 

y.  MRDLL.  alpha  is  the  multipath  delay  to  be  estimated.  C  is  the  noise 

y.  covariance  matrix . 

global  R  Beta  CM; 

y,  STEP  1. 

y,f=V_a.  This  is  the  function  that  we  want  to  be  minimized. 

H=sigl_H (Beta, alpha, M) ; 

f=(inv(C)*R) ’*(C-H*inv(H’»inv(C)*H)*H’)*(inv(C)*R) ; 

g(l)=-alpha; 

g(2)=alpha-1.6; 
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